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ABSTRACT 


An experimental study was conducted to investigate the use of 
advanced composite materials in unstiffened thin-walled cylindrical 
structures subjected to hydrostatic pressure. Specifically, graphite/epoxy 
was looked at for potential application for small submersible pressure 
hulls. The investigation concentrated on the fabrication, instrumentation, 
testing, and evaluation of several 12-ply graphite/epoxy cylinders of four 
different laminate configurations. 

The test specimens were 18 inches long with radii of 2.76 inches, and 
each was instrumented with 8 circumferentially and 2 axially mounted strain 
gages. Data collected from the strain gages throughout the tests were used 
in a modified Southwell type analysis and an analysis to deduce the 
buckling mode shape. The Southwell analysis method was shown to have a 
valid applicability to the buckling of laminated fiber-reinforced composite 
pressure vessels. A method was developed to use the bending strain data at 
locations around the centerline of the specimen to deduce the number of 
circumferential buckling waves. 

In addition, 1/2 inch ring specimen were cut from the excess ends of 
each specimen and tested under two-point compression, yielding information 
about the material properties of the laminates as fabricated. 

The hydrostatic tests were designed to verify the results of an 
analytical and numerical study that was conducted to determine the optimal 
laminate stacking sequences for buckling considerations. The test results 
indicate that the classical solution method presented by R.M. Jones 
produced accurate predictions of buckling strengths and correct relative 
rankings of the different laminate types. 
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NOMENCLATURE 


The following is a list of the special symbols and nomenclature used in 


this report. 


N 


Nye Ny 


Ai; 
€ 


K 


: 


= force resultants per unit length [1lbs/in] 
laminate extensional stiffness matrix [lbs/in] 
laminate extension/bending coupling matrix [lbs] 
laminate bending stiffness matrix [lbs in] 

= Strain {in/in] 


= Gurvature [1/in] 


My» My> My = moment resultants per unit length [lbs] 


u, 


<¢ mm ct WwW rms qQ 


Bes ce xz SD 


y, z 
v, W 
CLT 
GRP 
GFRP 
W/D 
AUV 
CTE 
FEA 


= stress [lbs/in-] 

= equivalent modulus [1bs/in‘] 

= mean radius [in] 

= wall thickness [in] 

= Young's Modulus [lbs/in’] 

= Poisson's ratio 

= hydrostatic pressure (positive external) [psi] 

= radial deflection (positive inwards) [in] 

= length of the cylinder [in] 

= number of circumferential buckling waves 

= number of axial buckling waves 
axial, circumferential and radial coordinates 
axial, circumferential and radial displacements 

= Classical Laminate Analysis 

= glass reinforced plastic 

= graphite fiber reinforced plastic (also CFRP) 

= weight to displacement ratio 

= autonomous underwater vehicle (also UUV) 

= coefficient of thermal expansion [in/in°F] 


= finite element analysis 
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1. INTRODUCTION 


1.1. Potential Benefits of Composite Pressure Hulls 


Since the 1960's fiber reinforced composite materials 
have been proposed for the construction of submersible 
pressure hulls. The potential benefits that can be offered by 
composite materials over metals include: 

o structural weight reductions; 

o fabrication of novel hull geometries; 

o incorporation of acoustic silencing; 

o electromagnetic signature reduction; 

Oo virtual elimination of corrosion. 
The purpose of this project is to gain insight into the use of 
graphite/epoxy for this application, through the systematic 
testing to hydrostatic pressure collapse and subsequent 
analysis of several thin-walled laminated composite cylinders, 
of various laminate types. 

Weight Reduction. Composite materials offer the 
Opportunity to make reductions in the portion of a 
submersible's weight that has to be devoted to the hull 
structure. In fact, submersible designers use as a measure of 
hull efficiency the ratio of the hull weight to the weight of 
the water the hull displaces. This ratio is referred to as the 
"W/D". It must be less than one for the submersible to float, 


and a smaller "W/D" means that the hull is more efficient, 
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Figure 1.1. Relationship of weight to depth ratio to collapse depth 
for several material systems. 


because for a given depth it can either carry more payload or 
use less energy. Values approaching 0.30 are required for the 
design of small vehicles with adequate allowances for mission 
capabilities. Figure 1.1 shows graphically the relationship of 
W/D to collapse depth for several materials [1]. From this 
figure it is seen that composite materials can offer 
substantial increases in collapse depth at any given W/D, or 
that for a specified collapse depth requirement composites 
enable large reductions in structural weight [2]. 

The major structural weight savings obtained via use of 
composite material construction can be directly utilized in 


other mission enhancing areas. Table 1.1 compares the typical 


weight allocations for an Autonomous Underwater Vehicle (AUV) 
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[1] with current U.S. Navy attack submarine (SSN) and fleet 
ballistic missile submarine (SSBN) design practices. Small 
submersibles are very volume critical, and designers always 
have more and larger payloads than they are capable to carry. 
Weight reductions in the hull structure can be translated into 
increased payload weight and volume capacity. For the many 
small submersibles where speed or range increases would be 
more beneficial, portions of the structural weight savings 
could be channeled into larger propulsion machinery or fuel 
Carrying capacity. For submersibles that require the 
capability to reach greater safe operating depths, the 
incorporation of composite materials could make the hull 
structure stronger without giving up established payload, 


propulsion and fuel allotments. 


30.6 
| 2. Variable Ballast | 3.7 | 12 | 12 | 











Table 1.1. Typical weight allocations of submersibles. 


Novel Hull Shapes. The methods of fabricating composite 
material hull structures make it easier for the designer to 


create complex geometrical shapes and incorporate advanced 
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Sandwich and double hull structural configurations. The use of 
novel structural geometries such as figure-eight and swedge 
hull shapes, allow for improved structural efficiency by 
optimizing the usage of space and volume as compared to 
traditional cylindrical and spherical hullforms [3]. 

Filament winding fabrication techniques allow for complex 
streamlined axisymmetric shapes to be constructed with 
relative ease. 

Signature Control. A major advantage of laminated 
composite pressure hulls is in the reduction of acoustic aes 
non-acoustic signatures. This can be crucial in missions such 
aS mine hunting and covert reconnaissance. Composites have 
seen years of usage as sonar windows due to their inherent 
acoustic transparency [4]. Also, by innovative utilization 
of low density sandwich cores and other acoustic damping and 
decoupling lamina and anechoic coatings, internal noises from 
machinery and electronics can be silenced to provide improved 
stealth. Because they are non-metallic, these materials are 
also non-magnetic and electrically non-conductive, greatly 
lowering a submersible's susceptibility to enemy weapons and 
sensors. 

Corrosion Sonia. Carencaion. which is a major problem 
for metals in the marine environment, can be virtually 
eliminated with properly selected composite material systems. 
However, in completeness, it must be mentioned that initial 


attempts at utilizing composites, specifically glass 


EZ 





reinforced plastic (GRP), found that seals around hull 
penetrations, hatches and viewports were susceptible under 
cyclic loading. The high deformation of the GRP as compared to 
metallic inserts left cut fibers vulnerable to exposure to 
water at high pressure. Because of this designs have been 
developed where cylindrical composite pressure hulls are used 
in conjunction with specially designed titanium end closures. 
In these designs all penetrations are made through the 
titanium end closures thus avoiding the difficult problem of 
penetrating the laminated composite pressure hull. Due to the 
weight of the metal end closures and because the unstiffened 
hulls are buckling critical, the overall hull weight may not 
show as significant reductions as predicted for all-composite 
hulls. Notwithstanding this, all the benefits presented by 
composite materials and the fact that any weight reduction is 
precious to a submersible designer, these advanced materials 
are well suited for utilization in pressure vessel design and 
construction. Graphite fiber reinforced plastics (GFRP) 
because of their higher modulus, lighter weight, and 
extraordinary fatigue properties generally offer more long 


term benefits than GRP. 
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1.2. Purpose of Hydrostatic Testing 


Four areas of major concern in the design of laminated 

composite pressure hulls are: 

o fabrication procedures, 

o repairability and serviceability, 

o inspection and certification, 

o in-service hull monitoring. 
The present investigation, consisting mainly of the 
hydrostatic testing of 12-ply graphite/epoxy cylindrical 
specimens of various layups, is aimed at providing data and 
techniques to better facilitate the addressing of these 
concerns. 

Prior to and concurrent with this test program, an 
analysis, both classical and numerical, was performed to 
predict the buckling characteristics of the graphite/epoxy 
cylinders that were tested here [5]. In this analysis, 
closed form solutions were able to be used for the cylinders 
constructed of orthotropic laminates. These laminates do not 
exhibit any coupling between extension and shear or bending 
and twisting stiffnesses (i.e. they have no Ax, Ar, Bis, Bre, 
Dis, Or Dog terms in their ABD matrices, which relate stress to 
strain). Finite element analysis was also conducted on these 
Orthotropic laminates and additionally on laminates with angle 
plies. These angle-ply laminates generally have some coupling 


between extension and shear or bending and twisting due to the 
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presence of these additional stiffness terms. 

Each specimen was instrumented with eight equally spaced 
circumferential strain gages. The average of these strains was 
compared to the axisymmetric strain that is predicted by 
theory and was also used as the membrane strain for a 
Southwell-type analyses to predict perfect shell buckling 
characteristics from experimental data from imperfect 
specimens. 

In summary, the main work on this thesis was to conduct a 
large series Bee oer nents in an effort to verify the results 
of theoretical elastic stability analyses methods that have 
been adapted for laminated cylindrical shells subjected to 
hydrostatic pressure. Since shell buckling characteristics are 
highly susceptible to laminate imperfections, geometrical 
deviations, and residual stresses, the results of this testing 
depended greatly on the effectiveness of the specimen 
manufacturing method. The fabrication technique chosen for 
this project consistently produced specimens of high 
circularity, but the following defects were found to varying 
extents in all the cylinders: wrinkled fibers, non-uniform 
resin distribution, voids and delaminations. In several of the 
specimens noticeable defects, especially severely wrinkled 
plies, were apparently responsible for large degradations in 
the buckling strengths of these cylinders. 

Although these circular cylindrical shells exhibit 


Significant buckling sensitivity to geometric shape 
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imperfections just like isotropic shells, in general, 
laminated composites offer the opportunity to reduce these 
effects by ser AEE Te appropriate fiber orientations. An 
overall goal of this entire analytical and experimental test 
program is to develop and verify a method to optimize the 
selection of ply orientations in order to achieve the most 
efficient hull for specific requirements. Before optimization 
routines can be applied to numerical or closed form solutions, 
there must be confidence that these solutions are fairly 
weeerate and can at least rank various laminate configurations 
correctly. Four different laminates were chosen to be 
evaluated by experimental test methods in an effort to test 
the analyses: 

o 1:1 specially orthotropic - [90;/0;]; , 

o 2:1 specially orthotropic - [90,/0,], , 

Oo quasi-isotropic - [90.,2/(+30).]s; , 

o interleaved (unsymmetric) orthotropic - [90,/0O]l«4, . 
By testing multiple specimens of each type, and comparing the 
rankings with the results of the analysis, insight into the 
correctness of the theoretical methods was obtained. 

Another significant contribution of this hydrostatic 
testing results from the extensive strain gage instrumentation 
that was utilized on each cylindrical specimen. On the 
centerline of every specimen were placed eight circumferential 
strain gages at 45° separation and two axial strain gages at 


180° separation. Readings from these gages provided the 


16 





ability to validate the assumptions of buckling prediction 
theories, to investigate the applicability of alternative 
methods of stability determination, and to provide insight 
into the mode shapes into which the cylinders buckled. 

In the 1930's, R.V. Southwell proposed a method to 
determine the buckling strength of an ideally straight column 
by analysis of test data taken from a column with initial 
curvature [6]. In the 1950's, G.D. Galletly and T.E. 

Reynolds proposed that a extension of this method may be used 
to determine the elastic stability of ring-stiffened cylinders 
[7]. The applicability of this method to determine the 
stability of composite laminated cylinders was investigated 
uSing the strain data collected during the hydrostatic 
pressure tests. 

By isolating the bending strain from the membrane strain 
at the eight locations around the centerline of the specimens 
and applying elementary concepts that relate the local radial 
deflections to the changes in circumference, a deduction was 
made of the buckling mode shape of each of the specimens. 
These mode shape deductions are compared to those assumed and 


predicted by the analysis methods. 
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2. ANALYTICAL BACKGROUND 


This section addresses the significant issues in material 
science and in structural Stability pertinent to the failure 
of laminated, fiber reinforced composite cylinders subjected 
to external hydrostatic pressure. Designing with laminated 
composite materials offers the opportunity to achieve specific 
structural properties through proper selection of: 

o the system of resin, fibers and enhancing agents; 

Oo processing and production techniques; 

Oo individual ply orientations. 

To make these selections the designer must consider concerns 
ranging from the micromechanical bonding of the fibers and the 
matrix, to the ability of the structural fabrication technique 
produce highly circular cylinders. This chapter attempts to 
bring together some of the information and experience that 
have been gained in these broad fields. Only by making the 
material scientist and the naval architect aware of the 
Other's knowledge and concerns, can all the benefits presented 
by the utilization of composite materials in submersible 
pressure hull design and construction be realized. 

The discussion first addresses the material science 
concerns, beginning with the micromechanical factors 
controlling the compressive failure of fiber reinforced 
laminated composites, and leading to a macromechanical 


treatment of Classical Laminate Theory. 
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The discussion then covers the major structural concerns, 
starting with the axisymmetric compression of circular 
cylinders, progressing through the compression of thin rings 
to the elastic stability theories for long cylindrical shells. 

The two fields merge with a coverage of a classical 
buckling analysis method for laminated cylinders, that is 
based on Donnell's stability theory. Figure 2.1 shows the 
coordinate system and dimensional nomenclature used in this 
analysis. Next, brief coverage is given of some recent finite 
element studies on this subject. 

Finally, a modified Southwell-type method is presented 
that should be universally applicable to a wide range of 
isotropic and anisotropic material systems. Additionally, a 
method is proposed whereby the bending strain data can be used 
to deduce the number of circumferential waves in the actual 


buckling mode shape of a cylindrical specimen. 





Figure 2.1. Diagram of the coordinate system used in the analysis. 
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2.1. Compressive Strength of Composite Laminates 


The compressive response of fiber-reinforced composite 
materials is still one of the least understood in the field of 
composites today. The compression failures of composite 
laminates can be categorized into the following mechanisms: 
"global Euler buckling, microbuckling, transverse tension, 
fiber kinking, fiber compression failures, matrix compression 
failures, and delamination." [8] The compression failure 
modes of fibers, matrices and unidirectional ee are 
illustrated in Figures 2.2 and 2.3, which were taken from 
reference ll. The predominant mechanism is determined by 
factors such as the individual properties of each of the 
constituent materials, the orientation of the plies within the 
laminate, the geometry of the specimen, the method of load 
application, any waviness and misalignment of the fibers, the 
presence of voids and discontinuities in the matrix, and any 
stress concentrations. 

If the fibers are weak enough in compression, they could 
fail due to compression before the onset of buckling and cause 
ultimate failure of the laminate. But unlike tension, the 
compressive Scmemcfen of most composites is limited more by the 
properties of the matrix than the fiber reinforcement. 
Microbuckling can begin to occur because the resin is soft, 
initiating ultimate failure, even when the fibers are not weak 


in compression. Other reasons the matrix may not be able to 
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provide adequate lateral support to the fiber are: proximity 
to a free surface, the formation of lateral splitting, a high 
void content, or environmental softening. Desirable matrix 
properties include a moderate amount of stiffness and 
toughness, but not an excess. Elastic microbuckling is the 
dominant failure mode when either the fibers have no initial 
misalignment or the matrix iS very soft. The low compressive 
strengths of fiber reinforced composites may also be the 
result of kink band formation, an occurrence that was 
predicted by theory and is now possible to observe 
experimentally. This type of compression failure results from 
fiber misalignment and involves plastic deformation of the 
resin. 

Microbuckling. Fiber buckling theories are based on the 
Stability of fibers ina flexible matrix. The fibers may 
buckle in either of two modes: (1) when the fiber volume 
fraction is greater than 30 percent, the fibers tend to buckle 
in phase (shear mode), and (2) when the fiber volume is less 
than 30 percent, the fibers tend to buckle out-of-phase 
(extensional mode) [9]. In a study of the effects of 
various constituent properties, fiber arrays, and bowed, 
fripended and misaligned fibers, Greszczuk concluded that the 
microbuckling models cannot predict the failure of non- 
idealized, actual composites. He stated further that low- 
modulus resins result in microbuckling failure, but when 


higher modulus resins are used the failure mode changes and 
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the compressive strength of the fiber reinforcement controls 
the failure [10]. 

Kink Band Formation. When the resin system is adequately 
stiff and there exists any fiber waviness, regardless of the 
mode of incipient failure, final failure occurs as the 
formation of kink bands. Chaudhuri conducted an analytical 
study to identify the basic failure mechanisms of 
graphite/epoxy cylinders with thick walls under hydrostatic 
compression [11]. He identified the dominant compressive 
failure mode for thick composite cylinders to be the 
"formation and propagation of fiber kink bands at the 
microscopic level, triggered by the fiber misalignment defects 
formed during the manufacturing process, leading to the shear 
crippling failure at the macroscopic level." He found that 
this kinking was the lowest energy compression failure mode at 
the microlevel of graphite/epoxy composites, and was limited 
by the initial fiber misalignment, the ultimate fiber strain 
and the two transverse shear moduli, G,, and G,, of the 
laminate. 

Chaudhuri's study is important because it shows that kink 
band formation is a distinct and lower energy failure mode 
than elastic microbuckling, not due solely to the effect of 
plasticity in the presence of initial fiber misalignment. 
Chaudhuri also noted that the ratio of hoop to axial layers 
did not strongly influence the compressive strength of the 


thick-walled (R/t<6) composite cylinders he investigated. 
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Figure 2.2. Compression failure modes of unidirectional composites 
according to Greszczuk [11]. 
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Figure 2.3. (a) Compression failure modes of fibers, matrices and 
unidirectional composites according to Hahn and Williams; (b) kink 
band geometry according to Chaudhuri [11]. 
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Higher toughness usually means fiber kinking will occur 
at lower strains, but lower toughness (i.e. higher modulus) 
usually means that there is less resistance to delamination 
propagation. Therefore, there must be a tradeoff between the 
delamination associated with higher modulus resins and the 
earlier (i.e. at lower strains) microbuckling initiation found 
with tougher resins. 

Shuart presented a linear analysis of short-wavelength 
modes shapes for off-axis and angle-ply graphite/epoxy 
laminates and a non-linear analysis of laminate stresses, out- 
of-plane displacements, and interlaminar shear strains 
[12]. He showed that significant interlaminar shearing 
Strains develop before short-wavelength buckling, and he 
described a compression failure criterion that included both 
these failure mechanisms. 

Even though fibers do not creep, all polymer matrices 
exhibit viscoelastic creep behavior which is affected by water 
absorption. Mechanical impact damage and exposure to high 
temperatures can aggravate this problem. Irion and Adams 
presented results that indicate that longitudinally reinforced 
graphite composites, in general, do not creep any significant 
amount. They found that transversely reinforced graphite 
composites exhibited long periods at high strains with no 
creep followed by abrupt jumps. Because few of their specimens 
attained steady-state creep and also due to the fact that the 


strain versus time plots displayed extensive amounts of 
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irregularities, they inferred that local failures were 
occurring almost continually, but after initiation the 
failures were quickly arrested [13]. 

Transverse tensile stress in unidirectional composites 
has been proposed aS a reason experimental compression 
failures occur at much lower loads then predicted by theory. 
Even though this stress is small, unidirectional composites 
have low transverse strength. 

In view of all this the expected failure sequence for 
laminated shells constructed from plies of unidirectional 
graphite fiber-reinforced composites with fiber volume 
fractions greater than 30%, and loaded in hydrostatic 
compression can be summarized as follows: "microbuckling ina 
shear mode (in-phase buckling) will occur, initiated by 
regions of local inhomogeneities (voids, stress 
concentrations, areas Of weak matrix, free boundaries), and 
lead to final failure in the form of fiber kinking or fiber 
overstressing in bending." [8] 

A major limiting factor on the compression strength of 
laminated composites is the technique of fabrication. 
Manufacture of the specimen has significant effect on the 
composite's compressive response because the presence of local 
inhomogeneities and defects, which are difficult to model and 
otherwise account for, influence the compressive failure 
mechanisms on the microstructural level more than any of the 


other stress state failure mechanisms. Therefore, the chosen 
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method of fabrication must result in a low void content (less 
than 2 %@), because otherwise the fibers located adjacent toa 
void will not have sufficient lateral support. Laminates must 
also be fabricated with very low amounts of fiber misalignment 
(compressive strength can be highly sensitive to fibers out of 
alignment as little as 2°). Some research suggests that 
processing with high fiber content (greater than 60%) may 
result in lower void contents and, therefore, more precise 
fiber alignment. 

The 5208 epoxy resin used in the fabrication of the 
cylinders for this present study has a high modulus, with 
failure strains of 1.4%. This low resin toughness has a 
detrimental effect on the progression of the failure. Studies 
have shown that compression failure of brittle laminates, such 
as T300/5208, is quite sudden, and that it is very difficult 
to arrest the resulting fiber kinking [14]. Sinclair and 
Chamis conducted compression tests on unidirectional T300/5208 
graphite/epoxy laminates and their results were divided into a 
high-strength group where failure was predicted by fiber 
compressive failure, a medium strength group that was 
predicted by flexure or delamination,and a low strength group 


that was predicted by Euler Buckling [15]. 
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2.2. Classical Laminate Theory 


Before discussing the elastic stability of composite 
pressure hulls, a brief overview of Classical Laminate Theory 
(CLT) is given. The concept of equivalent moduli for 
orthotropic cylinders is also developed to enable 
relationships designed for isotropic rings and cylindrical 
shells to be adapted, without excessive complexity, to the 
case of the laminated composite rings and cylindrical 
specimens tested in this research. 

For laminated orthotropic materials the constitutive 
relationship between stress and strain is complicated by the 
directional nature of the individual plies in the laminate. 
The elastic modulus for isotropic materials is replaced by an 
ABD matrix that contains information on the extension, shear 
and bending behavior of the laminate and the coupling between 
the extension, shear and bending. The normal presentation of 


this constitutive relation is 
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In general, the B,;, terms couple the bending and stretching 
effects. If the laminate is symmetric about the middle plane, 
there is no such coupling and thus B,,;=0. In specially 
orthotropic laminates, which are cross-ply laminates composed 
of plies oriented at 90° relative to each other, there is no 
coupling between bending and twisting effects, and, therefore, 
A,¢=A2,.=D,,.=D.,=0. 

To enable comparison of the relationships in the 
following sections for the axisymmetric compression and 
elastic stability of long composite cylinders with the 
experimental results, equivalent properties for an orthotropic 
cylindrical shell were developed. 

The membrane constitutive relations for an orthotropic 
material, where the bending response and extension/shear 
coupling are neglected, are 


28 





Ny Ai, Aiz 9 ex 
Ny? = |Ai2 A22 0 €y wh 
a O O Ag, a 


This is inverted to solve for the strains, 


such that 
ex 41, 4,2 O | IN, 
€y( = 1412 422 0 14N, oe 
O 0 0 a 
xy 66 XY 
where 
A22 
as, = 
Ay, Ag2 ~ Ai2 
~A,2 
aiz = 
Ay, Ao2 Ay2 (7) 
Ay; 
Qo2 = 


For the case of hydrostatic loading, 
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and thus 


iL 
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Ay, A22 ~ Ai2 


and the equivalent modulus in the hoop direction is therefore, 
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Similarly, the equivalent modulus in the axial direction is 
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and the equivalent Poisson's ratio is 
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These expressions can be used later with the experimental 
results, to check the spraCeivé properties of the cylindrical 
specimens as fabricated. Note that these relations only apply 
to the case of hydrostatic loading, for other loading 


conditions other equivalent moduli must be developed. 
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2.3. Axisymmetric Compression of Cylindrical Shells 


It has been said that axisymmetric yield failure is the 
most important mode of failure of externally loaded pressure 
vessels because a “Sensible design should eliminate failure 
due to instability." [16] This is lofty goal for thin- 
skinned pressure vessels because initial geometric 
imperfections on the order of a tenth of the thickness can 
cause a catastrophic loss of buckling resistance of pressure 
vessels when subjected to external water pressure. 

One of the earliest solutions for the axisymmetric 
deformation of circular cylinders under external loading was 
published by von Sanden and Gunther in 1920 [17]. Their 
Original work unfortunately had two signs interchanged in the 
formula for circumferential stress, which can give incorrect 
values if used unknowingly. Also, these authors did not take 
into account the beam-column effect and its associated non- 
linear deformations, that was later introduced by Salerno and 
Pulos [18] and refined by Ross [23] into the form 
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Ross gives the complete solution, but here is presented only 


the solution for our case of a circular shell element, 
stiffened by equal-strength frames located at its ends 


(x=0,L), and deformed symmetrically about the midspan, 


w=A, cosh] F, (x- ) |eos|F, (4-5) + A, sinh 





where 
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and A, and A, are functions of F, and F, that reflect the 
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characteristics of the end frame and the boundary conditions 


assumption that dw/dx=0 where the edge connects to the end 


frame (x=0,L). The variation of the hoop and axial stress 


along the shell is due to the radial contraction, the 


Poisson's effect of the axial compression, and the beam-column 
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effect. These are expressed in the following relationships for 
hoop stress and axial stress, respectively, where the last 
term is the beam-column effect and is positive in the inner 


fibers and negative in the outer fibers. 





Zz 
HoopStress: 0, = _Ew_ pRv ,_ Et {| wy, d‘w (17) 
K Chee (1 —-V eer Bb 
2 
AxialStress: o, = - Pk + 5 FE (2 , aoe (18) 
2c 2(1-v*) \ R* axe 


Ss 





2.4. Elastic Stability of Long Cylindrical Shells 


The discussion of the elastic stability of laminated 
composite shells begins with a review of classical buckling 
theories for isotropic shells of revolution. This paper 
addresses only cylindrical shells whose thickness t is small 
compared to the other dimensions of the shell, especially the 
radius of curvature. The discussion begins with a thin ring 
under symmetric two-point compression and works toward the 
theories for a cylinder of length L and radius R. 

The theory for thin isotropic circular rings under 
compressive loading was presented by S.P. Timoshenko in his 
classic text on elastic stability [19]. A summary of that 
analysis is given, and extended for the case of laminated 
composite rings. The relationships obtained here will be used 
in the analysis of the results obtained from the two-point 
compression testing of 1/2 inch ring specimens that were cut 
from the excess ends of each cylinder. 

Consider a isotropic circular ring of unit width, 
subjected to a two-point compression load P as shown in figure 
2.3. The bending moment at A (and at B) is denoted M, and the 


moment at an arbitrary location @ is 


M = M, + = (1-cos@) (19) 


The differential equation for the radial deflection of this 
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Figure 2.4. A ring of radius R compressed by two forces P acting 


along a diameter. 


ring at any point 6 due to this loading is 
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From the symmetry of the problem where 
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the general solution of equation 20 is found to be 
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Timoshenko then applied Castigliano's theorem to the strain 
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energy expression to get 


My = =* (2-1) (23) 


Of interest to us is the radial deflection at the point 


of load application (@=n/2), which from equation 22 is 
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where w < Q indicates the shortening of the diameter or 
deflection inward toward the center. 

For the case of a thin (R>>t) laminated composite ring, a 
two-point, symmetric compressive force has a negligible direct 
effect on the stress on the inner and outer surfaces at point 
A, and this stress, therefore, is almost entirely due to the 
bending moment resulting from that applied point force. 
Additionally the compressing ring is assumed to undergo 
cylindrical bending ( x = O ) and not anticlastic bending. 


Under these conditions, equation 1 can be simplified to 
M, . 


With cylindrical bending and negligible Poisson's effect, we 


Dy, D2 
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get that 
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M, = Dy)Xy (26) 


From the definition of curvature as the second derivative of 


the displacement we see that effectively, 
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Substituting this effective bending stiffness of the laminated 
ring, D..°eb = EI , into Timoshenko's formula for an isotropic 
ring, equation 24, and solving for D,. gives 


R?{P 


On the inner and outer surfaces the strain is related to the 


curvature, such that 


= eo 2 
SRG Is ae (29) 
where €,° is negligible for this loading, as explained above, 
and Zipper = t/2 on the surfaces. Combining equations 23, 26 
and 29 gives a relationship for the hoop strain on these 


surfaces in terms of the applied compressive load P, such that 
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Solving for D,, gives 


Dee +0.0908 46 (2) (31) 


Therefore, we now have a way to approximate the D,, of the ring 
with the slope of the load versus displacement plot, equation 
28, and with the slope of the load versus strain plot, 
equation 31. These experimentally determined D,,'s can be 
compared to the one obtained from CLT and the properties of 
T300/5208. The complete ABD matrices for each of the laminate 
types were developed in this way and are contained in Appendix 
E. This comparison provides insight into the actual bending 
stiffness properties of each of the laminate types as 


fabricated. 


Next, consider now how the radial deflection in the ring 
under the two-point, symmetric loading would be changed if a 
hydrostatic pressure of magnitude p was applied. This uniform 
external pressure produces a compressive force S = pR in the 
ring of unit width. This force combined with the symmetric, 
two-point loading considered above produces the deflection 


expressed by 
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This means that each term in w increases as the compressive 


force S increases, becoming indefinite when 
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In this case the critical compressive force occurs when n=2, 


therefore 
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This thin ring analysis can now be extended to the 
buckling of long circular cylinders under uniform external 
pressure. The critical pressure p,.,i:1 on an elemental ring of 
unit width of the cylinder is obtained by considering 
poisson's effects that were not a part of the isolated ring 


above, i.e. substitute E/(l-v’) for E and t°/l2 for I. 


Perit = 


ca a 


This calculated critical pressure is applicable as long as the 
compressive stress does not exceed the material's proportional 


iimit. 


One of the most significant contributions to this subject 


was the general theory of shell instability developed by W. 
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Koiter in 1945 [20}. In his thesis, Koiter demonstrated 
conclusively that initial geometric imperfections play a 
dominant role in reducing the buckling load of shell 
structures. This helped explain the scatter in experimental 
buckling test results and the general inability of researchers 
to experimentally attain theoretically predicted critical 
pressures. The greatest detriment to achieving theoretical 
buckling pressure predictions is imperfections in the shape of 
the shell. Koiter showed that in isotropic shells an 
perfection amplitude of just 10% of the wall thickness of 
the shell can result in as large as a 40% reduction in its 
buckling strength. This can be uncontrollable in thin-skinned 
shells with R/t>>10 like those studied here. More recent work 
by Koiter and others has led to theories where the effect of 
small but finite imperfections are included in the stability 
analysis. These give the designer the ability to account for 
geometric imperfections that can be assumed to be present in 
the as-manufactured pressure vessel. Noteworthy of these is an 
equation derived by Kendrick [21] which describes the 

stress in an isotropic cylinder due to eccentricities, such 


that 
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where C, is the amplitude of the initial deviation in a mode 
shape cos(n@). The most important term in this equation is 
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a ean (37) 


which shows that the hoop stress is directly proportional to 
the amplitude of the initial eccentricities in the cylinder. 
It also verifies that a failure envelope exists for a cylinder 
under external hydrostatic pressure because the stress becomes 


Singular as the critical buckling pressure is approached. 
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2.5. Classical Buckling Analysis of Composite Shells 


This section reviews some of the more prominent theories 
on the buckling stability of laminated composite cylindrical 
shells. These theories are referred to as "Classical" because 
they all relate the concepts of Classical Laminate Theory 
(CLT) to one or more of the classic elastic stability 
theories. The analysis begins with circular cylinders 
constructed of orthotropic layers and is based on work 
presented by R.M. Jones. Next, the analysis considers the 
introduction of angle plies and is based on work by Booten and 
Tennyson. Lastly, the analysis proposes the application of 
optimization techniques to laminate orientation selection and 
is based on the recent work of several different authors, most 
notably Sun. Please note that subscripts following a comma 
represent partial differentiation by that subscript. For 
example, the term N,, is equivalent to the expression QN,/dx. 

In 1968, Robert M. Jones presented an exact solution for 
the buckling of a circular cylindrical shell with multiple 
orthotropic layers under any combination of axial compression 
and lateral pressure [22]. Jones showed that previous 
attempts to predict the buckling behavior of Peciiaverad 
shells were erratic because they neglected the coupling 
between bending and extension, which his theory considered. He 
based his work on Donnell's stability differential equations 


for circular cylindrical shells, with an implied membrane 
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prebuckled shape and simply supported end conditions. 

By applying the Kirchhoff-Love hypothesis, the three- 
dimensional orthotropic elasticity relations were reduced to 
two-dimensional stress-strain relations for the k‘® layer. The 


variations in stresses in the k“ layer during buckling are 


d0% = On de, + Ox dé, 


k k k (38) 
do, = O12 de, 2 QO>> dé, 
k k 
Oty = O33 SY. 
and the variations in strains during buckling are 
S€, = €,+2K, 
(39) 


S€, = €,+ZK, 


BY xy = €3+ ZK, 


where €,, €,, and €,; are the variations in reference surface 
strains, and x:, x2, and x; are the variations in reference 


surface curvature, such that 


€,-u ,t+v (40) 


o - --w>, x,=-20 (41) 


Thus the variations in stresses in the k‘ layer can be 


expressed 
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See Cale, 2K.) 
These variations of stresses over the shell layers during 
buckling are then integrated to give the variations of forces 


and moments during buckling; 


ON, = Ay,€, + A,,€, + By,k, + ByK, oS 
ON, = A2€, + A,2.€, + BK, + Bok, (44) 
Deane. + Blak, (45) 

5M, = B,,€, + B,,€, + D,,*, + Dx, (46) 
5M, = 156. 26.0 Dek. + Do (47) 
(48) 


5M, = B,,€, + D33X; 


These force and moment resultants, which account for the 
coupling between bending and extension, only require four 
independent material properties, E,,*, E,y*, Gy", and viy, per 
layer. 


L.H. Donnell proposed the following stability 
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differential equations for circular shells subjected to 


combined axial compression and lateral pressure [23], 


QiN,. + ON, , = 9 (49) 


6N,y,, + SN, , = 0 


cl = _ 
— 8M xe + 8M yy, acy — OM yx, xy ~ OM yy + BON, + NAW x + NyW yy = 0 


Substituting into these equations, the variations of the 
in-plane forces per unit length and Eve variations of the 
moments per unit length, equations 43-48, together with the 
previous stress-strain and strain-displacement relationships, 
equations 40-42, gives three stability differential equations 
in terms of the displacements, u, v and w: 

et gw = 0 


an Sie Vane We 0 (51) 


L,,u+L,,v+L,,w = 0 


where L,;; are appropriate partial differential operators, see 


reference 22. 


There are sixteen possible force and geometric boundary 
conditions of which a combination of four constitutes a 
defined problem. Our main interest is in simply supported 


ends, such that, 


(52) 
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These four conditions are satisfied by the following buckling 


displacements: 


vs 15 cos| 
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ye 10, sin| sin (53) 


cos| 


“i i oR 
x _ Ix x 
“le Ea 


“~—ee” 
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where the overline indicates the amplitude. Placing these 
expressions into equation 5l, leaves the following 


relationship, whose determinant must be zeroto obtaino be 


nontrivial solution: 


Py Fi. Fi 


c=] 


Po1 P52 P53 


mm\2 --/ n\2 
P35, F532 {Fs 2) -W{3) } it 


vt = 0 (54) 


This results in the following stability criterion, which we 


shall refer to as "Jones' Equation": 


W,( 2) +H, (3) =- FP +P Fy3F 2 7-F ii F 23 ae Digs Big (55) 
Pio y\p S3a 123 2 S é 
Teese: ae Ia IE. Fie 


where 
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(61) 








Anz , 2822 ay + =a (ay 
R L 


The buckling load under hydrostatic pressure is obtained by 


Substituting the following load factors: 


ee re 
N,= 25 = 2 (62) 


Ct 
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Jones' contribution was instrumental because it 
introduced the effect of the coupling between bending and 
extension that exists due to the presence of the multiple 
orthotropic layers in laminated, composite shells. The more 
difficult questions of initial geometrical and material 
imperfections and prebuckling deformations were not 
considered. 

In 1978, M. Booten and R.C. Tennyson conducted a study 
that considered the behavior of laminated anisotropic 
cylinders with small axisymmetric shape imperfections 
[24]. Their analysis used exact boundary conditions and 
nonlinear prebuckling deformations, and comparisons were made 
with experimental data. These authors concluded that the 
effect of varying the fiber angle on the interactive buckling 
behavior of laminated cylinders with "perfect" geometry and 
"less than perfect" axisymmetry is significant. In a later 
work with Hansen [25], Tennyson demonstrated that the 
buckling strength of an axially loaded cylindrical shell could 
be increased significantly through a judicious choice of 
laminate configuration. This has naturally lead to the topic 
of optimization. 

Composite materials offer the designer the ability to 
tailor material properties to the specific application. In 
filamentary composites, such as graphite-epoxy, this involves 
ply angle selection. where designers aspire for the capability 


to select optimal ply angle distributions for specific 
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geometries and load conditions. In the past ten years several 
studies have been conducted into laminate optimization for 
cylindrical shells. 

In 1981, Y.S. Nshanian of the then U.S.S.R. and M. Pappas 
examined the optimization of the ply angle distribution across 
the thickness of composite shells for buckling and vibration 
considerations [26]. They proposed methods to determine 
the optimal ply angle variation through the thickness of 
symmetric angle-ply shells of uniform thickness. USing a 
mathematical programming algorithm, the procedure attempted to 
maximize the buckling load of a thin, sSimply-Ssupported, 
circular, cylindrical, angle-ply shell. Their results showed 
significant gains using optimal variable ply angle 
configurations, aS compared to ones with an optimal constant 
ply angle. But reference shells with variable but non-optimal 
angle-ply configurations (0°/90° and 90°/0°) were often near 
optimal in performance and always superior to optimal 
conventional angle-ply shells. 

In 1987, Sun and Hansen attempted to maximize the 
buckling load of a laminated-composite, circular-cylindrical 
shell as a function of ply orientation [27]. They used a 
two-step optimization routine with a first step austen search 
to determine an initial guess and a second step systematic 
search based on work by Powell [28]. The optimization 
routines were applied to the finite difference solutions to 


von Karman-Donnell non-linear shell strain-displacement 
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relations. Experiments were conducted on selected, optimal, 
four-ply, graphite/epoxy cylinders and compared to certain 
reference laminates. The experimentally determined buckling 
load under external pressure for the optimized cylinder was 
3.19 times as high as the [0/90}, reference cylinder, but no 
[90/0}, reference cylinders, with hoop plies on the outer and 
inner surfaces, were tested under external pressure. 

In 1989, Sun presented another optimization study using 
linear theory and the same two-step optimization method 
{29}. The study assumed the cylinders were constructed of 
N orthotropic plies with the same thickness, with the i‘ ply 
containing an equal number of fibers in the +@, and -@,; 
directions uniformly distributed through the ply's thickness. 
This eliminates the stiffness coefficients of the laminate Ai,, 
A2>, Bis, Boe, Dig, and D2... This assumption is valid for 
laminates with many +@ and -@ plies of thicknesses much 
smaller than the thickness of the laminate. It serves only as 
an approximate result for laminates with relatively few plies, 
such as the 4-ply ones Sun tested and the 12-ply laminates 


tested here. 
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2.6. Finite Elements Analysis of Composite Shells 


In 1984, Abu-Farsakh and Lusher investigated the buckling 
behavior of glass-reinforced cylindrical shells under combined 
axial compression and external pressure [30]. Initial 
imperfections, load-deflection response, and buckling modes 
and loads were determined experimentally for 20 specimens and 
compared to the results of a finite element analysis (FEA). 
The FEA predicted lower buckling loads than were 
experimentally measured. They concluded that the size of the 
mesh must produce a higher number of divisions 
circumferentially than the number of circumferential buckling 
waves to obtain reasonable results. 

As a prelude to the hydrostatic testing outlined in this 
report, extensive finite element analysis work was completed 
by T. Perry at the University of Massachusetts, Lowell, 
Massachusetts [5]. The suitability of three finite element 
software packages were evaluated for use in this project, EMRC 
NISA, Algor's SUPERSAP, and HKS's ABAQUS. ABAQUS, which has 
anisotropic material capability for buckling and non-linear 
analysis, proved to be the best of the three packages. ABAQUS 
finite element models for specially orthotropic, orthotropic, 
quaSi-isotropic, and isotropic cylinders were run using the 
material properties of T300/5208 graphite/epoxy and simply 
Supported boundary conditions. The models were created using 


16 circumferential by 10 longitudinal (16cx101) parabolic 5 
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degree of freedom/node elements for a total of 2560 degrees of 
freedom (dof). The FEA gave stiffer results by as much as 25% 
over the classical analysis uSing Jones' equation described 
above. Trends were evident but they were sometimes contrary to 
those predicted by Jones' analysis. The present experimental 
test program aims to determine which method, classical 
analysis or FEA, produces the best results, both in accurately 
predicting buckling strengths and in the relative ranking of 


different laminate types. 
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2.7. Southwell-type Analysis of Elastic Stability 


In 1932, R.V. Southwell, in an attempt to find "the 
practical value of a theory of elasticity," proposed a 
procedure to analyze test data from an elastic column with 
initial curvature in an effort to determine the buckling 
characteristics of a perfectly straight column [6]. This 
procedure was an attempt to validate theories of elastic 
stability experimentally, realizing that "(l) the finite 
strength of actual era ecmaanid (2) the unavoidable 
imperfections of workmanship" prevent the practical 
realization of elastic stability's concept of a "critical 
load." In the method, the central deflection of the column was 
plotted against a load-normalized deflection and the slope of 
the best straight line drawn through the test data was 
interpreted as the critical buckling load of the a perfectly 
straight column. Southwell foresaw that the basic concepts of 
this method should be widely applicable to structural 
Stability problems wherever "inaccuracies in the specimen or 
in the experimental apparatus introduce displacements which 
increase continuously with the load." 

Recognizing its wide applicability to elastic stability 
problems, in 1955, Galletly and Reynolds developed an 
extension of Southwell's method that could be used to 
determine the elastic general instability pressure of ring- 


stiffened cylindrical shells subjected to external hydrostatic 
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pressure [7]. In this procedure the applied hydrostatic 
pressure was plotted against the circumferential bending 
strain that existed largely due to the initial out-of- 
roundness of the test specimens. 

L.H. Donnell concluded back in 1938 that Southwell's 
straight-line plot method was probably not applicable to 
curved structures because the deformed shape of curved 
structures was not a developable surface and the corresponding 
differential equations would be nonlinear [31]. Galletly 
and Reynolds showed that although the plot of pressure versus 
the reciprocal of the bending strain due to out of roundness 
was not linear, it was so slightly curved as the buckling 
pressure was approached that it could be closely approximated 
by a Straight line. Another paper, by Horton and Cundari, 
demonstrated using previously published experimental data and 
theoretical analysis that the method was "of absolute 
generality," and they concluded that "theories for perfect 
bodies ... not in agreement with experimental values derived 
on the basis of Southwell type plots must be considered 
inadequate, incomplete or in error." [32] 

Galletly and Reynold's method was based on small- 
deflection theory where displacements and strains grow 
indefinitely large as the buckling pressure is approached. 
Because strains are easier to measure experimentally than 
displacements and intercepts are easier to objectively 


determine than asymptotic slopes, they plotted the reciprocals 
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of the strain quantities versus pressure rather than 
displacement versus pressure like Southwell had originally. 
An Pe oieanat ion of the method of applying Southwell's 
concepts that was used to analyze the strain data from this 
test program is presented graphically in the next three 
figures. Figure 2.5 shows some typical circumferential strain 
data plotted versus the applied hydrostatic pressure. The 
curved lines show the variation of the circumferential strain 
with pressure at three individual points on the circumference 
Be a cylindrical specimen (0°, 45°, and 90° relative to each 
other). These curves become hyperbolic as the critical load is 
approached, but notice that they can be approximated by 
straight lines for a large portion of their length, as shown 
by the dashed lines. The slopes of these straight lines are 
different at each different location around the cylinder due 
to the varying amounts of local eccentricities. Galletly and 
Reynolds illustrated that the average of a large number of 
evenly distributed circumferential strain gage readings 
closely approximates the axisymmetric circumferential strain 
predicted by the Equation 13. This membrane strain can be 
considered as the elastic circumferential strain that would 
exist in a perfect structure whose radius is being uniformly 
compressed. The individual circumferential strains that 
actually occur on the cylinder are due to a combination of 
this membrane strain and the bending strain induced by the 


initial local eccentricities in the circular cylindrical shape 
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Figure 2.5. Typical plot of circumferential strain versus applied 
hydrostatic pressure, showing extraction of bending and nonlinear 
strains. 


of the structure. We therefore consider the average of the 
experimental strains as the membrane strain, Enxeatrane, and the 
difference between this and the actual individual strain gage 
readings as the local bending strain, Eensing, aS Shown in 
Figure 2.5. 

A plot of the bending strain versus the applied pressure, 
p, iS nonlinear and Evening ASYMptotically approaches infinity 
at certain values of p, the first of which is the minimum 
elastic buckling pressure, Pp... The other values of p where 
E€,ending DECOMEeS infinite do not occur in reality and therefore, 
a plot of the reciprocal of the experimentally measured 
bending strain, 1/€sening, VerTSusS the applied pressure will be 
curve that intercepts the p-axis at pui:, the minimum buckling 
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Bieure 2.6. Plot of 1/€ponding versuS pressure to determine a 


cylinder's critical buckling pressure. 


pressure of a perfect cylinder. 


Referring again to Figure 2.5, the gages located at the 


relative locations of O° and 45° behave in a manner that makes 


this method of determining p.ui. feasible. However, the fact 
that the strains from the gage located at 90° cross the 
membrane strain line prevents the utilization of data from 


this gage in this exact method. This is because where € equals 


Enembrane the value Of 1/€,..14::, Decomes infinite. The Southwell 


plot of all three gages where Ea,average = Encabrane 1S ShOwn in 


Figure 2.6. 


One way to be able to use the readings from all "good" 
gages is to plot the reciprocal of only the nonlinear 


component of the local bending strain, instead of the bending 
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mreure 2.7. Plot of 1/€non-linear versus pressure to determine a 
cylinder's critical buckling pressure. 


strain itself. In Figure 2.5, the strain at each gage is 
noticed to be very linear for the first 50 to 75 psi of 
pressure. The dashed straight lines are extensions of this 
linear portion of the strain readings. We can, therefore, 
break each bending strain response curve into a linear portion 


and a nonlinear portion we will call Enoniinesar, SUCH that: 


7 (63) 
©hending = CP * €noniinear 


Where c is the constant slope of the dashed line. A plot of 
€.,onlinear VETSus the pressure also, as expected, becomes infinite 
at the critical buckling pressure. The three gages of Figure 


2.5 were used in this manner to create the plot of 1/Enoniinear 
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versus p in Figure 2.7. 

When one recognizes that the above deduction makes no 
assumptions or references to the material characteristics of 
the cylinder, the applicability to the difficult problem of 
laminated fiber reinforced composite cylinders is realized. 
These methods were used in an attempt to validate previous 
theoretical and finite element analysis of laminated composite 
cylinders, using strain data collected from the hydrostatic 
test of composite cylindrical specimens. Predicting the 
critical buckling pressure during a test before the onset of 
failure is also a possibility with this method. It is proposed 
that a method could be developed whereby strain data froma 
pressure hull loaded to only a percentage of its predicted 
buckling load could be used to predict how that cylinder will 
behave at pressures approaching the critical buckling 
pressure. Development of a method of this type would enable 
for hull certification and in-service monitoring to determine 


degradations from impact damage, cyclic loading, etc. 


39 





2.8. Buckling Mode Shape Analysis 


A method is proposed to determine the buckling mode shape 
of a circular cylinder using the baSic concepts that R.V. 
Southwell presented for relating the radial displacement of a 
point on a axisymmetrically compressed cylindrical tube to its 
change in circumference [33]. For this analysis, the 
radial displacement is assumed to occur in a plane 
perpendicular to the longitudinal axis, and the circular 
cross-section is assumed to remain ates, The fractional 
contraction of the circumference results in a new radius (R- 
w), and can be expressed as an axisymmetric circumferential 


strain, 


beeciCh-w)>2mk | (64) 


y 2uR 


by | x 


This can be converted into an expression for the radial 
displacement in terms of just the nominal radius R and the 


uniform circumferential strain €,, 


w= -€,R (65) 


The actual circumferential strain at any point on a cylinder 
being compressed nonuniformly can be broken into a membrane 
stain that would define the uniform contraction in the radius 
and a bending strain that would be the local contraction or 
extension relative to this uniform change. Therefore, by 
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looking at just the affect of the bending strain at several 
points around the circumference, an idea of the number of 
circumferential buckling waves can be deduced. 

This method is applied to each of the specimens using the 
eight circumferential strain gage readings at the last 
pressure where all gages were recorded prior to buckling 
failure. To better illustrate the deformed shaped the radial 
deflections were magnified from 200 to 500 times normal size. 

These mode shape plots are contained in Appendix C. The 
odd-numbered figures RoR Pain the variations of bending strain 
with pressure at the noted locations around the circumference 
of each specimen and the even-numbered figures contain polar 
plots of the deduced mode shapes with the magnified radial 
deflections. 

Application of the above analysis to these bending 
strains and magnifying the radial displacements, even though 
it does not give absolute deflections, does make the buckling 
mode shape clearly apparent. This experimentally derived mode 
shape can then be used directly in the classical and finite 
element buckling analyses, or as a verification of the 


analyses' circumferential buckling wave assumptions. 
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3. EXPERIMENTAL PROCEDURE 


The main contribution of this thesis is the results 
obtained from the hydrostatic pressure testing to failure of 
cylindrical test specimens of various laminate configurations. 
The purpose of this testing was to experimentally verify the 
results obtained from classical and finite element buckling 
analyses. In addition to checking how precisely the analysis 
methods predict the buckling failure pressures of these 
miitively Simple laminates, these methods must be verified to 
at least rank the buckling pressures of the different laminate 
types correctly. The correct ranking 1s essential to selecting 
one of these analysis methods as a basis for studies into the 
optimization of ply orientations for buckling strength. 

The test matrix originally called for the testing of at 
least three 12-ply composite cylinders from each of the four 
laminate configurations. In all, thirteen cylinders were 
tested but because of fabrication induced defects only eight 
specimens were considered "good." The additional 5 defective 
specimens were tested to failure, but their strain data and 
comparative value were Regee Leda anese cylinders had 
noticeable flaws and would normally have been discarded prior 
to testing, but in order to refine the test set-up and also in 
an effort to gain an understanding of how these defects affect 
the failure mechanisms, all were tested, but only the strain 


data from the good specimens was used in the analysis and 
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presented in the results. 

The strain gage readings taken throughout the hydrostatic 
tests and the critical buckling loads measured comprise the 
bulk of the experimental data. The results were tabulated and 
plotted in various formats, and then analyzed to determine 
patterns or conclusions about the strength of composite 
cylinders subjected to external hydrostatic pressure. The 
experimental results were also compared with the analytical 
results to check predictions about laminated cylinder buckling 
strengths and assumptions of buckling mode shapes. 

Very little research has been conducted to experimentally 
determine the behavior of thin-walled composite cylinders of 
various lay-ups subjected to hydrostatic external pressure to 
failure, and to compare the results with theoretical 
predictions. Currently there is only a meager amount of test 
data on the buckling of thin-walled composite shells under 
hydrostatic pressure. The data collected and analyzed here is 
hoped to give a better understanding of the ability to use 
composites for submersible hull structures. Ultimately, 
insight is sought into improved ways to select material 
combinations and laminate configurations that will enable the 
benefits of laminated composite material technology to be 
maximized in future design and construction of submersible 


pressure hulls. 
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3.1. Cylindrical Test Specimens 


All of the tested specimens were fabricated from 12-ply 
laminated circular cylinders with nominal diameters of 5.52 
inches, skin thicknesses of 0.072, and lengths of 18 inches. 

All the cylinders were constructed of individual plies of 
T300/5208 graphite/epoxy pre-impregnated unidirectional ply 
(prepreg uniply). Table 3.1 contains the nominal properties of 
T300/5208 in this format and the typical properties of various 
other types of graphite fiber / polymer matrix unidirectional 
composites, for comparison. These properties were obtained 
from the Composite Materials Research Group at the University 
of Wyoming and are for nominal 60% fiber volume composites. 
Digital micrographs of polished specimens cut from the excess 
ends of the laminated tubes tested here indicated an average 


fiber volume of 50%. 
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Table 3.1. Typical Material Properties. 


A. Fabrication Method 

There are several different viable fabrication methods 
that could have been used to produce the cylindrical specimens 
for the hydrostatic testing. The primary methods used in the 
construction of tubes and cylindrical vessels out of laminated 
composite materials are hand layup with vacuum bag curing, 
filament winding, and tape winding. Hand layup with vacuum bag 
curing was not considered to Ne eaecuate because it would 
almost be impossible to control wrinkles and geometric 
deviations from perfect circularity. As discussed earlier, 
Slight imperfections have a significant effect on the buckling 


failure characteristics of cylindrical shells in compression. 
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Figure 3.1. The elastomeric tooling fixture used to manufacture the 
cylindrical test specimens. 


These irregularities could be largely eliminated and at least 
controlled by constructing the specimens with filament winding 
technique. This fabrication method produces highly circular 
and wrinkle-free cylinders, but distribution of the fibers 
could not have been made as uniformly as required. Tape 
winding, which combines the best attributes of hand layup and 
filament winding was a technically feasible option, but the 
cost of a tape winder was prohibitive to project budget 
constraints. For pneae reasons, a novel method of constructing 
the cylindrical specimens with high circularity and minimal 
laminate defects had to be incorporated. An elastomeric 
tooling method that involves the use of an expanding internal 


mandrel and a highly circular external mandrel was chosen to 
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fabricate the test specimens. 

This method calls for an internal mandrel constructed of 
a high coefficient of thermal expansion (CTE) material which 
expands in radius during the increased temperature curing 
cycle. This mandrel is wrapped with thin sheets of elastomer 
to give a consistent, smooth surface to the final cylinder. 
Sheets of prepreg uniply are then hand laid up onto the 
mandrel. The sheets of uniply must be carefully measured, cut 
and rolled onto the cylinder by one person while the another 
person makes sure that the lamina lays seat smoothly and at 
the desired angle. After all plies are applied to the internal 
mandrel, two highly accurately machined semi-circular metal 
external tools are clamped on. The pressure exerted during the 
curing process by the expanding high CTE internal mandrel 
forces the multiple layers of uniply tape up against the 
Stiff, low CTE external tool. 

This method proved to consistently create laminated 
cylindrical specimens with extremely close circular 
tolerances. The details of the elastomeric tooling fixture 
used to manufacture the specimens for this research is 
diagrammed in Figure 3.1. It has an internal mandrel made of a 
Beatactic foam core (CTE = 28 x 10° in/in°’F) that is wrapped 
with several layers of silicone rubber sheet (CTE = 162 x 10° 
in/in’F), and an external tool fabricated from two 1/2" thick, 
24" long aluminum pipe halves (CTE = 14 x 10° in/in°F) which 


were bored to a .005" circular tolerance. 
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B. Laminates Tested 
The cylinders each had twelve plies and were constructed 
out of one of the four following different layups: 
[90,/0,], - 2:1 Specially orthotropic 
[90;/0,;], - 1:1 specially orthotropic 
[90,/(+30).],; - quasi-isotropic 
{[90,/0]., - interleaved orthotropic 
The numbers indicate the orientation of the individual plies, 
the numerical subscripts indicate the number of plies of that 
orientation, the alphabetic subscript S indicates that the 
laminate is symmetrical about the mid-surface, and the 4T 
indicates that the unsymmetric laminate is constructed of 4 
identical plies. Each of these laminates are diagrammed in 


Figures 2.2, 2.3, 2.4, and 2.5. 


68 





90 
90 


PPP 


i 
4 


[90,/9,]. 


alll 1:1 Specially 
0 


90 


0 


| 


2:1 Specially 


i. Orthotropic L Orthotropic 
ee. 90 
90 %0 
% a 
[a rz 
peeeure 3.2. Bigune. suo. 
90 
| 90 
ll 
[ 90,/ @30), ]. ~([90,/0] ,. 


2:1 Interleaved 
Orthotropic 


Quasi-isotropic 





rj 
fate 
as 
a} 
a 
ee) 
Nn 


Bigure 3.4. 





The 1:1 and 2:1 specially orthotropic laminates do not 
exhibit any stiffness coupling and, therefore, are the easiest 
of the anisotropic laminates to analyze. This also means that 
the material does not warp or twist during post-cure cooling. 
This is not as much of a problem with cylinders as flat 
plates, but was a concern nonetheless. 

The unsymmetric, interleaved orthotropic layup, which is 
commonly used in industry, has somewhat better interlaminar 
properties than the specially orthotropic laminates because 
the perpendicular ply interfaces are distributed through the 
thickness. This causes a reduction of the interlaminar 
stresses that cause delamination, a major failure mechanism in 
laminated composites. 

The quasi-isotropic layup was used to compare the 
composite cylinder failure characteristics to cylinders made 
from isotropic materials. Laminates made with this type system 
of ply orientations result in directional stiffness 


Characteristics very close to isotropic. 
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3.2. Hydrostatic Pressure Tests 


A. Test Set-up and Instrumentation 

The hydrostatic testing was conducted in an Atlas 600 psi 
air pressure chamber located in the Mechanical-Thermal Lab at 
C.S. Draper Laboratory in Cambridge, MA. Pressure was 
provided to the chamber via a pressure control system that was 
accurate to within +2 psi. The cylinders were potted at both 
ends in specially designed endcaps. The endcaps, shown in 
figure 2.6, were designed with a slot for potting that has an 
angled outside edge so that the edges remained circular, and 
also so that only the outside edge of the tube would actually 
be in contact with the endcap. Armstrong C7 epoxy was used to 
pot the tubes in the endcaps. This epoxy was compliant enough 
to allow the ends of the cylinders to rotate when loaded, 
Simulating simply supported edge conditions. 

Each specimen was instrumented along the centerline with 
10 Micro Measurements 4 inch, 120 2 strain gages, model number 
EA-06-500BH-120. Eight of the gages were evenly distributed at 
45° separations around the centerline of the specimen, and 
oriented to read circumferential or hoop strain. The other two 
gages were positioned on the centerline opposite each Bier at 
a 180° separation, and oriented to read axial strain. The gage 
leads were draw out of the pressure chamber, via breakout 
boxes, to a ten channel BLH Model 1225 switching and balancing 


unit, and then to a BLH Model 1200 digital strain indicator. 
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Figure 3.6. Details of cylinder potted in endcap. 


The test set-up is diagrammed in Figure 3.7 and a photograph 


of the facilities is provided in Appendix E. 


B. Test Procedure 

The testing of each cylindrical specimen was conducted by 
the following procedure: 
o The specimen was placed in the chamber, the strain gage 
leads were connected to the breakout box, the gages tested for 
integrity, and a copper bleed-off line was connected to one of 
the endcaps. 
o The pressure chamber was closed and sealed, and each of 
gages was manually zeroed with the switching and balancing 
unit. 
o The pressure in the chamber was increased to 25 psi and 
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Figure 3.7. Diagram of experimental set-up, showing pressure chamber 
and instrumentation. 


maintained long enough for strain readings to stabilize and to 
be manually recorded for each individual gage. 
o The pressure was increased to 50, 75, 100, 110, and so on at 


10 or 5 psi increments until failure, taking readings at each 


step. 
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3.3. Ring Tests 


The external tools that were used in fabrication of the 
cylinders were 24 inches long but the specimens were only to 
be 18 inches long. This left 3 inches of excess on each end of 
the cylinder. Several 1/2 inch rings were cut from the excess 
end of each cylinder, and these were used to help determine 
some of the in situ material properties of the specimens. 

Two strain gages were applied at the same angular 
location of each ring, one on the inside and ays other on the 
outside surface, 90° from the point of load application. The 
rings were then loaded in two point compression on the 
Instron® 4505 located in the Composite Materials Laboratory of 
the Material Science and Engineering Department at the 
Massachusetts Institute of Technology, Cambridge, MA. 
Measurements of strain and displacement versus load were 
recorded throughout each of the tests. The data acquisition 
system consisted of a four channel Micro Measurements strain 
gage conditioner unit, model 1220A, and a National Instruments 
Lab-PC integrated data acquisition card. 

These ring tests gave some additional mechanical property 
data for the various laminate configurations in the "as 


fabricated" condition. 
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4. RESULTS and DISCUSSION 


4.1. Ring Tests 


Throughout the two-point, symmetric compression test of 
the each 1/2 inch ring, load, strain and displacement data was 
collected. Plots were made of the load versus displacement and 
the load versus strain data. The slopes of the load versus 
displacement curves were used with the relationship of 
equation 28 to calculate a "w test" D,, for each of the ring 
specimens. Additionally, the slopes of the strain versus 
displacement curves were used with the relationship of 
equation 31 to calculate an "€ test" D,.. These experimentally 
determined D.,'s were compared with the D,, from CLT contained 
in the ABD matrices of Appendix E, which were developed using 
the material properties for T300/5208 given in Table 3.1. This 
comparison is summarized in Table 4.2 below. The experimental 
"w-test" and "“€-test" values are seen to agree well with each 
other, and are in reasonable agreement with the values 
calculated using CLT. 

Because the radius, thickness and width of the ring 
specimen has a significant impact on equation 31, these 
dimensions were accurately measured with a micrometer at eight 
locations around the ring and the average of these values used 
in calculations. These averaged values for R, t and b are 


presented in Table 4.1. 


es: 















6] 2.759 | 0.072 | 0.523 | 
[saver | 0.070 | 0.503 | 
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| Specimen R 15 b 
No. fin. ] Cairns | fin. 4 | 










w-test | €-test CLT 
| 10100_| 361.6 | 352.1 | 526.2 | 
| 10520_| 367.5 | 377.4 | 574.4 | 


Table 4.2. Predicted and Measured D>. 





134.9 
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4.2. Cylinder Tests 


A total of 13 full-length cylindrical specimens were 
tested to failure, but several of these had noticeable defects 
and thus only served to refine the test procedure and gain 
more experience with the graphite/epoxy pressure vessels with 
large defects. Table 4.3 contains a breakdown of the specimens 
tested and their critical buckling pressures. This table also 
gives comments on each "bad" specimen, explaining the details 
of Eric problem, such as any noticeable flaws. In all, eight 
specimens have been tested that were classified as "good". 
Only the buckling pressures and strain data from these eight 


specimens were analyzed for this investigation. 
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The original endcaps were designed with large cones to 
seat the cylinder edges onto and were constructed out of 
steel. Reproducible end conditions were extremely difficult to 
create using these steel endcaps. Specimens 5 and 15 were the 
first two cylinders tested and the conical endcaps were used. 
The cones on the endcaps seemed to wedge into the cylinders 
under pressure creating erroneous results. These problems 
inspired the development of a new endcap design. The details 
of the potting of the cylinders in these endcaps are shown in 
Figure 3.6. Because of the change in endcaps the results from 
the first two cylinders tested were discarded. 

Specimen 10 was a very "dry" laminate and had excessive 
leakage through the plies themselves. This specimen was coated 
in C-7 epoxy to prevent the pressure from equalizing inside 
the cylinder and enable testing to be completed. Because of 
this coating, the results were not considered accurate. 

During the manufacture of specimens 7 and 8 the pre-cured 
laminate that was rolled onto the internal mandrel was too 
large for the aluminum external tools to be clamped around. 
These specimens were put into a freezer to induce shrinking 
and then forcibly clamped into the tools. Both these specimens 
had creases of the order of the laminate thickness from the 
pinching of the semicircular aluminum tools. Failure of these 
Specimens under hydrostatic loading was premature and both 
fractured cleanly down the creases as shown in the photographs 


in Appendix F. For this reason these specimens were not used 
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for comparative purposes. 

Table 4.4 summarizes the critical buckling pressures that 
were obtained from Jones‘ equation, ABAQUS modelling, 
experimental testing, and Southwell-type analysis for the 
"good" specimens. 

Although the test matrix was not as large as originally 
planned, "good" data was obtained for two 1:1 specially 
orthotropic, two 2:1 specially orthotropic, three quasi- 
isotropic, and one interleaved (unsymmetric) orthotropic 
specimen. Because of the close similarity between the data 
from each of the "good" specimens from the same layup, these 
results can be confidently used to evaluate the analysis 


methods. 
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Table 4.4. Predicted and Measured Critical Buckling Pressures. 


A. Critical Buckling Pressures 


As stated earlier, a major objective of this work is to 
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verify which of the analysis methods not only best predicts 
the buckling characteristics, but consistently correctly ranks 
the stability of each of the laminates used. This latter 
objective is crucial as these methods are proposed for further 
studies into optimizing the orientations of the plies ina 
multi-layered composite cylinder to achieve the highest 
buckling strength. 

Figures 4.1, 4.2, 4.3 and 4.4, contain bar graphs which 
each present the ranking of the four different laminate types. 
Figure 4.1 gives the values obtained using Jones' buckling 
criterion, equation 56, which was derived from Donnell's 
Stability equations and CLT. Figure 4.2 are the values given 
by the FEA program ABAQUS. Figure 4.3 presents the average 
critical pressure for each laminate type predicted by the 
Southwell-type analysis, taken from the plots in Appendix B. 
Figure 4.4 gives the average experimental buckling pressures 
for each laminate type as obtained from the hydrostatic tests 
of the "good" specimens. 

These figures show that the ABAQUS finite element 
analysis program does not rank the laminates correctly. This 
is a major obstacle to using ABAQUS to optimize the selection 
of individual ply orientations. The ABAQUS analysis predicts 
that the orthotropic laminates with two hoop plies for every 
axial ply will have the greatest buckling strength. It 
predicts that the 2:1 symmetric specially orthotropic laminate 


1s strongest and that the interleaved, unsymmetric orthotropic 
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laminate will have a higher critical buckling pressure than 
either the 1:1 specially orthotropic of the quasi-isotropic 
laminates. This was not collaborated by the testing. 

The Jones' equation analysis predicted that the l1:l 
specially orthotropic laminate would give the cylinder the 
greatest buckling strength, followed by the 2:1 specially 
orthotropic, the quaSi-isotropic, and then the interleaved 
orthotropic laminates, respectively. This ranking was exactly 
the same order as the critical pressures measured in 
hydrostatic tests and predicted by the subsequent Souchwell- 
type analysis. This is a significant finding. Applying 
optimization techniques to the classical analysis methods of 
Jones will be much less expensive than to apply similar 
techniques to the finite element methods of the ABAQUS 
program. That these classical methods have been shown to 
correctly rank at least the popular laminate types tested 
here, means that designers can apply them to optimization 


studies with a substantial amount of confidence. 
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Figure 4.1. Laminate ranking using Jones' Equation. 
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Figure 4.2. Laminate ranking using ABAQUS. 
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Southwell Analysis Predictions 
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Figure 4.3. Laminate ranking using Southwell-type 
analysis. 
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Figure 4.4. Laminate ranking using the results of 
hydrostatic testing. 
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Next, the actual critical pressures from the hydrostatic 
tests are compared to the buckling pressure predictions to 
verify the accuracy of the various analysis methods. Figures 
4.5, 4.6, 4.7 and 4.8 contain bar graphs which each present 
one laminate type and show a quantitative comparison of the 
different buckling pressures predicted by the three analysis 
methods and obtained from hydrostatic testing. 

The ABAQUS program was very accurate in its prediction 
for the cylinders constructed out of the 1:1 specially 
orthotropic and the quasi-isotropic laminates, but the 
magnitude of the overestimation of the buckling strength for 
the 2:1 specially orthotropic and interleaved, unsymmetric 
orthotropic laminate cylinders prevents this finite element 
method from being used with any confidence to predict the 
critical pressures of angle ply laminated cylinders. 

Jones' equation gives predictions of the buckling 
strength larger than those measured experimentally, but within 
reasonable deviation. It is significant to note the precision 
and accuracy with which the Jones' equation and the average 
Southwell-type analysis exhibit, relative to each other. Table 
4.5 contains the results of these two analyses which both 
consider the strength of a "perfect" geometry circular 
cylinder. It should be noted that the Jones' predictions were 
for 60% fiber volume laminates, and recently taken digital 
micrographs of the specimens indicate an average 50% fiber 


volume. Recalculations using the updated laminate properties 
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show that the values obtained from the two methods are 


actually in even better agreement. 


Type Equation Plot Difference 
| interleaved orthotropic | 164 | _155_—| + 5.8 8 


Table 4.5. Comparison of Jones' Equation and Southwell Plot. 
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Figure 4.5. Critical buckling pressures for 1:1 Specially 
Orthotropic laminates. 
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Figure 4.6. Critical buckling pressures for 2:1 Specially 
Orthotropic laminates. 
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Figure 4.8. Critical buckling pressures for Interleaved 
Orthotropic laminates. 
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B. Axial Strain Response 

It was hoped that the two axially mounted strain gages 
could be used to investigate the beam-column effect of the 
long cylindrical specimens. The plots of the axial strain 
responses for all of the good specimens are contained in the 
even-numbered figures of Appendix A. The axially strains for 
the 1:1 specially orthotropic specimens seemed to be quite 
linear with little or no increase in curvature as the critical 
pressure was approached. The plots of the axial strain 
response for the 2:1 specially orthotropic, quasi-isotropic, 
and interleaved (unsymmetric) orthotropic specimens showed the 


more traditional hyperbolic shape at increased pressures. 


C. Circumferential Strain Response 

Plots of the circumferential strain responses from the 
Strain gages located on each of the specimens are contained in 
the odd-numbered figures of Appendix A. All the cylinders had 
eight circumferentially-mounted strain gages evenly 
distributed around the centerline at 45° separations. Note 
that when one of the plots has less than eight circumferential 
strain curves, it means that one or more of the gages failed 
during the test. ; 

The hyperbolic curvature of the circumferential strain 
response measured by each of the gages, as predicted by 


theory, is clearly evident in all cases. The individual 


circumferential strain responses to increasing hydrostatic 
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pressure can be broken into a membrane and a bending strain 
response. The bending strain results from the shape change of 
the circular cylinders' cross-section as the buckling 
pressure is approached. 

This excellent circumferential strain data is a 
significant contribution of this research, and proved well- 
worth the time-intensive labor it required to apply the gages, 
and then to collect and analyze the extensive amount of data 
that was yielded. Because the eight circumferential Strain 
gages were evenly distributed around each cylinder's 
centerline, they provided the necessary data to construct 
detailed Southwell-type plots and to perform an analysis to 


deduce the buckling mode shape for each cylinder. 


D. Average Circumferential Strain 

By averaging the readings of the eight circumferential 
gages at each pressure we get a close approximation of the 
axisymmetric circumferential membrane strain of each of the 
cylinders as they were compressed. This is the strain response 
that would exist if the cylinders cross-section would remain 


circular throughout the loading. 
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Figure 4.9. Axisymmetric circumferential membrane strain 
response for the 1:1 specially orthotropic cylinders. 
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Figure 4.10. Axisymmetric circumferential membrane strain 
response of the 2:1 specially orthotropic cylinders. 
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Figure 4.11. Axisymmetric circumferential membrane strain 
response for the quasi-isotropic cylinders. 
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Figure 4.12. Axisymmetric circumferential membrane strain 
response for the interleaved orthotropic cylinder. 
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To examine how accurately the averages of the eight 
circumferential strain gage readings approximate the membrane 
strain, the experimental average of the circumferential strain 
gage readings is compared to the theoretical axisymmetric 
circumferential membrane strain response predicted by equation 
9 for hydrostatic loading. The curves for each specimen and 
the theoretical curve for each laminate type, are contained in 
Figures 4.9, 4.10, 4.11 and 4.12. 

These plots show that the average of the eight 
Mrecierential strain gage readings from specimens made from 
the same laminate, are almost exactly the same. These average 
strains also compare extremely well with the axisymmetric, 
Circumferential membrane strain calculated using the material 
properties for T300/5208 with 60% fiber volume and equation 9, 
which is based on CLT. The slight discrepancy could be the 
result of the fact that the laminates actually had only a 50% 
fiber volume, not 60% as assumed in the material properties. 
Additionally, equation 9 does not consider the bending 


response or the extension/shear coupling in the laminate. 


E. Modified Southwell Analysis 

Using the method described in the section on Southwell 
analysis, where the reciprocal of the non-linear strain has 
plotted versus pressure, Southwell-type plots were constructed 
for each good specimen. The odd-numbered figures of Appendix B 


show the Southwell plots for each of the specimen over the 
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same large range of pressure. The even-numbered figures of 
Appendix B show the details of the Southwell method by 
fleplaying only a narrow range of concern for each specimen. 
The dotted lines on these plots are the extensions of the 
reciprocals of the non-linear strain curves. Close examination 
of these plots show the extraordinary small amount of scatter 
that was involved in the prediction of critical buckling 
pressures. The values of these predicted buckling pressures 
for the specimens had they been of perfect geometry is given 
in the tables and bar charts earlier in this chapter. As noted 
before, the predictions of Jones' equation were consistently 
only between 5 to 7.5 percent higher than those obtained from 


the Southwell plots. 


F. Mode Shape Deduction from Bending Strains 

Using the method outlined in the section on Mode Shape 
Analysis, the bending strains of each specimen were utilized 
to deduce the buckling mode shape of each cylinder. The plots 
of the bending strains and the subsequent variations in shape 
are contained in Appendix C. The buckling mode shapes are 
drawn at the last pressure where all the strain gages were 


read and recorded. 
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Table 4.6. Predicted and Experimentally Deduced Mode Shapes. 
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It is significant that each and every specimen, 
regardless of laminate type, appeared to buckle ina 
distinctive n= 3, three lobe mode shape. This was predicted 
by Jones' equation for all cases, whereas the FEA predicted 
n=2 to be the lowest energy buckling mode shape for both the 


1:1 and 2:1 specially orthotropic cases. 


G. Failure Progression 

The 5208 epoxy matrix used in the laminates tested in 
this investigation is a very brittle resin, having a low 
toughness and high modulus (1.4% failure strain). Brittle 
resins tend to be unable to resist the propagation of 
delamination, exhibiting quite sudden and uncontrollable 
failures. 


The failure of these T300/5208 laminated cylinders was 
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exceptionally catastrophic, appearing to happen in 
microseconds. Attempts were made to film the collapse of the 
cylinders with a high-speed (250 frames per second) video 
camera, but even when slowed to one frame per second for 
viewing the failure sequence was too fast to capture on film. 
All the good specimens shattered spectacularly, leaving 
only very small pieces of intact laminate, large quantities of 
individual broken fibers, and basically only a dust residue of 
the resin. Both the specimens that were creased during 
fabrication (cylinders 7 and 8) fractured cleanly along the 
crease, leaving large sections of the cylinder relatively 


intact. 
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5. SUMMARY AND CONCLUSIONS 


The hydrostatic tests conducted for this thesis have been 
successful in achieving the original objectives. 

The classical and finite element buckling analysis 
methods were examined for accuracy and correctness in 
Getermining the critical buckling pressure of unstiffened 
laminated composite cylinders of various layups. From this 
examination, Jones' equation for the buckling of generally 
orthotropic cylinders has proven accurate in predicting the 
buckling loads and correct in the ranking of different 
laminate types. The validity of this equation which has been 
experimentally verified, and its relative simplicity make it 
well-suited for use in optimization programs that could be run 
economically on personal computers. 

A systematic test program was developed whereby 
cylindrical specimens can be hydrostatically tested to 
buckling failure under repeatable conditions. The 
instrumentation used yielded valuable strain data that was 
used in Southwell-type analysis and also in determining the 
buckling mode shape of the specimens. 

Comparison oe the Southwell plots and the Jones' equation 
predictions, showed the consistency with which this type 
analysis of circumferential strain data from the test of 
geometrically imperfect specimen can be used to predict the 


buckling strength of perfect cylinders. The applicability, 
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which had previously been shown for isotropic ring-stiffened 
cylinders, has proven to be applicable to unstiffened 
cylinders constructed of laminated anisotropic materials as 
well. 

Through careful documentation of the fabrication process 
and through the utilization of a large diversity of material 
characterization tests, insight into the role fabrication 
induced imperfections, residual and interlaminar stresses play 
in the degradation of pressure hull performance has been 
gained. In addition, ring compression see have shown that 
the bending stiffnesses predicted by classical laminate theory 


may for symmetric laminates be exaggerated. 
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6. RECOMMENDATIONS 


This testing and analysis has shown the large increase in 
knowledge that can be gained by the extensive strain gage 
instrumentation. It is recommended that in future testing the 
pressure be steadily increased and the strains be recorded 
continuously on a computer data acquisition system. In these 
tests, the pump on the chamber would kick on and off with 
quite violent jolts attempting to maintain a steady pressure. 
This disrupted the peapiication of the strain gages and 
caused delays before readings could be recorded. The effect of 
these shocks and delays on the tubes, which is unquantifiable, 
would be eliminated with continuously increasing pressure 
loading. 

It is believed that by modifying the Southwell and mode 
Shape analysis procedures presented here that a practical 
method could be developed whereby the compressive buckling 
strength of individual pressure hulls could be qualified. A 
properly instrumented pressure hull, tested to only a 
percentage of its predicted collapse load, would exhibit 
strains and displacements that may enable the non-destructive 
verification of predicted collapse loads. 

This test series should be extended to angle-ply 
laminates to check Jones' equation's accuracy for non- 
orthotropic ply orientations. Notwithstanding the lack of 


angle-ply testing conducted here, it is recommended that 
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Jones’ equation be utilized in laminate optimization 
investigations in lieu of finite element methods, which have 
shown to improperly ant the buckling strength of even the 
Simple laminates tested here. 

Graphite/epoxy laminated thin cylindrical shells such as 
those tested here, have shown to have potential in the design 
of small submersibles that are only required to operate at 
shallow depths. It is recommended that in order to avoid 
catastrophic failure of the hull structure, and subsequent 
complete i of the vehicle, that utilization of tougher 


resins be investigated. 
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APPENDIX A 


CIRCUMFERENTIAL AND AXIAL STRAINS 
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Figure A.1. Circumferential Strains around the centerline 
of cylinder 6. 
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Figure A.2. Axial Strains at the center of cylinder 6. 
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Figure A.3. Circumferential Strains around the centerline 


of cylinder 9. 
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Figure A.4. Axial Strains at the center of cylinder 9. 
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Figure A.5. Circumferential Strains around the centerline 
of cylinder 12. 
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Figure A.6. Axial Strains at the center of cylinder 12. 
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Figure A.7. Circumferential Strains around the centerline 
of cylinder 13. 
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Figure A.8. Axial Strains at the center of cylinder 13. 
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Figure A.9. Circumferential Strains around the centerline 
©f cylinder 16. 
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Figure A.10. Axial Strains at the centerline of cylinder 
16. 
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Figure A.11. Circumferential Strains around the centerline 
of cylinder 17. 
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Figure A.12. Axial Strains at the center of cylinder 17. 
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Figure A.13. Circumferential Strains around the centerline 
of cylinder 18. 
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Figure A.14. Axial Strains at the centerline of cylinder 
18. 
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Figure A.15. Circumferential Strains around the centerline 
of cylinder 23. 
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Figure A.16. Axial Strains at the center of cylinder 23. 
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APPENDIX B 


SOUTHWELL-TYPE ANALYSIS OF STRAINS 
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Figure B.1. Southwell Plot for Cylinder 6. 
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Figure B.2. Detail of Southwell Analysis for Cylinder 6. 
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Figure B.3. Southwell Plot for Cylinder 9. 
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Figure B.4. Detail of Southwell Analysis for Cylinder 9. 
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Figure B.5. Southwell Pliot for Cylinder 12. 
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Figure B.6. Detail of Southwell Analysis for Cylinder 12. 
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Figure B.7. Southwell Plot for Cylinder 13. 
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Figure B.8. Detail of Southwell Analysis for Cylinder 13. 
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Figure B.9. Southwell Plot for Cylinder 16. 
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Figure B.10. Detail of Southwell Analysis for Cylinder 16. 
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Figure B.11. Southwell Plot for Cylinder 17. 


0.007 
0.005 
O7803 
07001 
=¢ 00) 


= 005 


| Ve e nonlinear 


=—O0> 





Power t g 
150 160 170 180 190 200 


Pressure (psi) 
Figure B.12. Detail of Southwell Analysis for Cylinder 17. 
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Figure B.13. Southwell Plot for Cylinder 18. 
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Figure B.14. Detail of Southwell Analysis for Cylinder 18. 
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Figure B.15. Southwell Plot for Cylinder 23. 
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Figure B.16. Detail of Southwell Analysis for Cylinder 23. 
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APPENDIX C 


BENDING STRAINS AND MODE SHAPE PLOTS 
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Figure C.1l.. Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 6. 
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Figure C.2.. Shape variation due to bending strain for 
Cylinder 6. 
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Figure C.3.. Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 9. 
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Figure C.4.. Shape variation due to bending strain for 
Cylinder 9. 
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Figure C.5.. Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 12. 
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Figure C.6.. Shape variation due to bending strain for 
Cylinder 12. 
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Figure C.7.. Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 13. 


Cylinder No. 13 
(90570-1. 90° 





Peay = 173 psi 


Zi 


*Radial deflections shown at 1 70 psi, magnified 200 times. 
Figure C.8.. Shape variation due to bending strain for 
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Figure C.9.. Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 16. 
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Figure C.10.. Shape variation due to bending strain for 
Cylinder 16. 
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Figure C.11.. Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 17. 
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Figure C.12.. Shape variation due to bending strain for 
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Figure C.13.. Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 18. 
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Figure C.14.. Shape variation due to bending strain for 
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Fiqure C.15.. Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 23. 
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APPENDIX D 


SPREADSHEETS OF STRAIN GAGE READINGS 


iL 43) 





STRAIN GAGE 


cylinder # 


6 


DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 


layup test date 


[903/03jJs 2/19/92 


1. Raw Circumferential Strain Gage Data. 


Pres Art 

0 8 
ZS -80 
50 -124 
75 +-136 
100 ~-160 
110 -14]1 
120 -147 
130 -148 
140 -148 
150 -151 
160 -149 
165 -141 
170 -141 
i795 -137 
180 -132 
183 


Be 
-2 
86 
29 

0 

-33 

-35 

-47 

-61 

-71 

-91 

-~96 

-95 

-105 
-110 
-118 


2. Intercept at 


Ar 
INT - 36 
3. Zeroed 

Pres Asx 
0 0 
25 -44 
50 -88 
75 -100 
100 -124 
110 -105 
320 -121)1 
130 -112 
140 -112 
150 -115 
160 -113 
165 -105 
170 -105 
175 -101 
180 -96 


183 


Be 
143 


C 

3 
-174 
-297 
-413 
= 339 
Sino 
(55) 
=655 
—1o3 
-775 
-680 
- 668 
=701 
Sie 
-849 


P=0 using 


C 
=o 


Circumferential 


Be 
0 
-57 
-114 
-143 
-176 
-178 
-190 
-204 
-214 
-234 
-239 
-238 
-248 
-253 
-261 


G 

0 
-123 
-246 
ee 
-488 
-524 
-580 
-644 
~702 
-724 
-629 
=—Ol7 
=650 
-704 
-798 


D E F G H 
-8 2 6 3 -11 


= 1S): -~145 =e eye -88 -144 
-~239 -257 =295 a pe: = 76 16 
= 367 ~343 -428 -170 -~314 


-445 -410 -568 -181 -387 
-468 -427 -634 -177 ~422 
oi -~439 -690 -161 -438 
=> 50 =455 -767 -145 - 466 
-576 ~465 -852 -118 -490 
-627 -472 -995 -73 =o20 
-645 -466 -1056 -41 =os0 
S16 5174 -452 -1104 -14 -542 
-672 -44) -1181 24 =593 
-687 -~427 -1243 70 =o) 
-711 ~401 -1338 115 -564 


values at 25 and 50 psi. 


D E F G H 
= 7. 35 =15 = 31 -52 


Strain Data (using intercept). 


D E F G H 
0 0 0 0 0 
-108 -~112 -140 -57 -92 
-216 ~224 -280 -114 -184 
-344 =3 10 -413 -139 -262 
- 422 -377 255 -150 -335 


-445 -394 =-619 -146 =370 
-482 ~406 -675 -130 - 386 


-527 -422 -752 -114 -414 
=553 - 432 -837 -87 -438 
-604 =4.39 -980 -42 ~475 
=622 -433 -1041 -10 -486 
-629 -419 -1089 17 -~490 
-649 -408 -1166 55 =501 
-~664 -394 -1228 101 =505 


-688 “368 -1323 146 S05) 72 


JPME 


] 
-189 
= Sis 
-510 
-675 
~784 
-829 
2) 

-1030 
=1 201 
=1275 
=—ls36 
=1425 
=1521 
216133 
=1730 


~26 


-163 
=326 
-484 
-649 
-758 
-803 
-899 
-1004 
= eS 
-1249 
-1310 
-1399 
-1495 
-1587 
-1704 


9 
=65 
-82 
-68 
3518) 


oF 

81 
35 
225 
276 
321 
380 
448 
a1? 


-54 


-14 
~28 
—14 

21 

56 

91 
i332 
189 
277 
330 
375 
434 
502 
571 





STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 


cylinder # layup test date 
9 {[903/03)s 10/31/91 
1. Raw Circumferential Strain Gage Data. 
Pres Ar Be C D E F G H 
0 -22 -17 -18 3 ] -15 -23 -18 
25 -104 -98 - 312 -145 -184 -263 -147 -~156 
50 -137 -149 -416 -259 -284 -472 -213 - 260 
75 -189 -197 -444 - 382 -~412 -673 -288 - 376 
100 -247 -256 -513 -454 -526 -856 - 339 -471 
110 -245 -275 -535 -470 -562 -936 -344 -506 


-1035 - 358 -547 


120 =-274 -293 =974 -496 -613 
-1134 = 359 =-583 


S30 -286 -316 -610 ~a 17 =656 


m0 -291 -333 -651 =932 =689 -1240 307 2020 
moO 8-308 -358 -714 = 212), -734 -1375 =-345 O59 
160 -318 -387 -788 -~554 -762 -1564 = s02 -708 
m0 =-327 -403 -868 =o30 -773 -1786 -230 -748 
meee 328 -422 -1017 -489 -757 -2077 -104 -832 
190 -325 -443 -1472 -265 -649 -2743 408 -1071 
200 -4279 2152 -2583 


2. Intercept at P=0 using values at 25 and 50 psi. 


As Bs Cc D E F G H 
INT -7)1 =47 -208 ~-31 -84 -54 -81 -52 


3. Zeroed Circumferential Strain Data (using intercept). 


Pres As Bs Cc D E F G H 
0 0 0 0 0 0 0 0 0 
ZO - 33 -51 -104 -114 -100 -209 -66 -104 


50 -66 -102 -208 —-226 -200 -418 -132 -208 
75-118 -150 -236 -351 - 328 =619 -207 - 324 
100 -176 -209 -305 - 423 -442 -802 =258 -419 
110 -174 -228 oe? - 439 -478 -882 -263 -454 
120 -203 -246 =366 -465 -529 -981 -277 -495 
130 -215 -269 -402 - 486 -572 -1080 -278 -531 
140 -220 -286 -~443 - 501 -605 -1186 -276 -564 
150 -237 -311 -506 -519 -650 -1321 -262 -607 


160 -247 -340 =200 -523 -678 ~-1510 -221 =6506 
170 -256 -356 -660 -507 -689 -1732 -149 -696 
180 -257 -375 -809 -458 -673 --2023 -23 -780 
190 -254 -396 -1264 -234 -565 -2689 489 -1019 
200 =422> 2235. =2531 


{2] 


J 

-21 
-194 
- 343 
- 496 
-645 
—699 
=-769 
-836 
Sel 
-1006 
= bs, 
~1234 
es 2 
N63 


~45 


-149 
-298 
=a I 
-600 
-654 
-724 
-791 
-866 
-961 
-1074 
- 1189 
- 1347 
-1608 


K 


a0 
SVE 
-324 
-381 
=390 
- 406 
-415 
-419 
-421 
a0 2 
- 350 
-251 

65 


-28 


-102 
-204 
2010 
=335 
-362 
-378 
- 387 
- 391 
=J05 
-374 
=e 
=223 

hes 


mn ane - 


-~ 


a 





STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 


cylinder # layup test date 
12 (904/02)s 11715791 


1. Raw Circumferential Strain Gage Data. 


Pres Ar Be Cc D E F G H 
0 -5 2 i7178 -14 -6 1 -4 ] 
25 -104 -94 17194 ~62 -~115 -~94 -96 -130 
50 -151 -142 17247 -66 -211 -137 -171 -395 
75 -203 -208 17298 -69 -354 -206 -269 -671 
100 -244 -247 17346 -58 - 459 -148 -303 -924 
moO. =-257 =%+-250 17370 -36 ~497 -151 -296 -1018 
120 -274 ~-260 17390 -18 -573 -170 -~298 -1188 
130 -289 -260 17414 10 -661 -187 -274 -1364 
140 -296 -247 17434 37 -733 -202 -242 -1538 
150 -306 -225 17453 81 -836 -214 -185 -1783 
160 -323 -174 17474 163 -1010 -251 -42 -2229 
165 -322 -123 17486 246 -1104 -264 66 -2475 
170 -332 -32 17500 389 -1334 -313 205-3012 
175 -349 54° 17511 543) =—1522 -377 530 -3522 
180 -419 2545 17521 985 -1889 -565 1182 -5200 
181 -5800 


2. Intercept at P=0 using values at 25 and 50 psi. 


Ax Be Cc D E F G H 
INT -57 -46 -58 -19 -51 -21 130 


3. Zeroed Circumferential Strain Data (using intercept). 


Pres Ax Be Cc D E F G H 
0 0 0 0 0 0 0 0 
25 -47 -48 b ~4 -96 ~43 ~75 -~265 
50 ~94 -96 a -8 -192 -86 -150 -~530 
75 -146 -162 d -11 ~335 -155 -248 -806 
100 -187 -201 0 -440 -97 -282 -1059 
110 -200 -204 g 22 -478 -100 -275 -1153 
120 -217 -214 a 40 -554 -119 -277 -1323 
H32Q@ -232 -214 g 68 -642 -136 -253 -1499 
140 -239 -20) e 95 -714 -151 -221 -1673 
150 -249 -179 139 -817 -163 -164 -1918 
160 -266 -128 221 -991 -200 -21 -2364 
165 -265 -77 304 -1085 -213 87 -2610 
170 -275 14 447 -1315 -262 306 -3147 
175 -292 100 601 -1503 -326 551 =-3657 
180 -362 297 1043 -1870 -514 1203 --5335 
181 -5935 


[42. 


-2 
-71 
-120 
-18] 
=202 
-198 
-198 
—LoZ 
-174 
-146 
= 


bo2 
at3 
840 


-22 


-49 
=98 
sles, 
-~180 
-176 
-176 
-170 
-152 
-124 
~49 
Pe, 
174 
335 
862 


~84 
- 168 
=276 
- 306 
- 306 
= Se 
=316 
=309 
he 
=253 
-208 
-123 
= 37 
210 


eg ee 





STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 


cylinder layup test date 
13 (904/02)s 10/17/91 


1. Raw Circumferential Strain Gage Data. 


Pres Ax Be Cc D E F G H 
0 -4 -6 1 -2 -5 -] -6 -5 
25 -120 -85 -114 -102 -122 10243 -85 -215 
50 -192 -139 -186 -202 -289 -182 -121 -435 
75 -267 -174 -249 -261 -366 1009 -107 -633 
100 -306 -163 -316 -309 -483 1131 -80 -962 
110 -314 =-160 -357 -328 ~542 1100 -53 -1123 
120 -325 -146 -397 -340 -~606 1133 -11 -1302 
130 -341 -123 -~451 -~341 -696 1133 48 -1549 
140 -364 -95 -516 -320 -836 1167 146 -1858 
150 -378 -30 -595 -271 -1024 1201 291 -2298 
160 -376 105 -~732 -138 -1311 1355 568 -3010 
170 -371 459 -1093 316 -1947 1602 1348 -4699 
173 -7500 


2. Intercept at P=0 using values at 25 and 50 psi. 


Ax Ba C D E F G H 
INT -48 -3)1 -42 -2 45 -49 5 


3. Zeroed Circumferential Strain Data (using intercept). 


Pres At Be Cc D E F G H 
0 0 0 0 0 0 0 0 
25 -72 -54 -72 -100 -167 b -36 -~220 
50 -144 -108 ~144 -200 -334 a ~72 ~440 
75 -219 -143 -207 -259 -411 d -58 -638 


100 -258 -132 -274 -307 -528 -31 -967 


MEO) -266 -129 =315 326 -587 g -4 -1128 
m0) -277 -115 =395 ~3368 -6351 a 38 -1307 
130 «2-293 ~92 -409 =~ 339 -741 g 97 -1554 
140 -316 -64 -474 -318 -881 e 95) —1863 
180 -330 1 =203 -269  -1069 340 -2303 
OO -328 136 -690 -136 -1356 617 ~-3015 
oO 6-323 490 ~-1051 318 -1992 1397 -4704 
173 -7505 


1343 


5 
-69 
-116 
-80 
-38 
a 
35 
94 
205 
367 
659 
1378 


eee 


15 


Os: 
Ze7 
389 
681 
1400 


4 
-109 
-~226 
-251 
-304 
—3 1s 
“315 
352 
-348 
= Ol 
-354 
-274 


-~117 
-234 
Pa), 
= Sa 
=320 
- 326 
-340 
=396 
- 369 
=362 
ri2oe 


ES 


——— 


ee eee oe A 





STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 


cylinder # layup test date 
16 (902/(+_30)2)s 11/5791 
1. Raw Circumferential Strain Gage Data. 
Pres A Fe B C D E G H 
0 } 4 10 -15343 6 10 13 -8 
ZO -23 -43 -136 -15504 -132 -139 -186 -105 
50 -98 -33 -303 -15574 -310 = 0 -385 -218 
75 -132 -1 -428 -15625 -445 -431 -610 -325 
100 -161 24 -546 -15676 -583 -564 -861 -391 
110 -158 58 -577 -15700 -618 -606 -976 -404 
120 -156 65 -589 -15718 -646 -635 -1097 - 391 
130 -164 27 -592 -15743 -671 -662 -1267 -376 
140 -166 -52 -561 -15779 -712 -675 -1541 -322 
150 -124 -124 -387 -15825 -691 -588 -1912 -154 
160 7 -269 269 -15915 -581 -227 -2674 350 
165 -5200 
2. Intercept at P=0 using values at 25 and 50 psi. 
Ae Fa B C D E G H 


31 5454 46 32 13 8 


intercept). 


INT a2 oS 


3. Zeroed Circumferential Strain Data (using 


Pres Ak Fe B Cc D E G H 
0 0 0 0 0 0 0 0 
25 -75 10 -167 b -178 -171 -199 -113 
50 -150 20 -334 a -356 -342 -398 -226 
75 -184 52 -459 re} -491 -463 -623 -333 
100 -213 Ty -577 -629 -596 -874 -399 
110 -210 1) -608 g -664 -638 -989 -412 
120 -208 118 -620 a -692 -667 -1110 -399 
130 -216 80 -623 g -717 -694 -1280 -384 
140 -218 ] -592 2 -758 -707 -1554 -330 
150 -176 -7) -418 -737 -620 -1925 -162 
160 -45 -216 238 -627 -259 -2687 342 
165 -5213 
170 
175 
180 
183 


34 


-184 
= 30 
=536 
-721 
-794 
-868 
=9o9 
=1095 
-1241 
-1449 


-177 
-354 
-529 
-714 
-787 
-861 
-952 

-1088 

-1234 

-1442 


-7 
=185 
-347 
-514 
-687 
-745 
-815 
=6972 

= 1027, 
-1187 
-1545 


=23 


-~162 
=324 
-491 
-664 
ee 
~io2 
-869 
- 1004 
-1164 
alles 77 





STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 


cylinder # layup test date 
17 (902/(+_30)2)s 11/8/91 

1. Raw Circumferential Strain Gage Data. 
Pres Ar Be C D E 2 G H J 

0 -49 -146 -42 -16 -15 8 -13 -47 = 3) 
25 -92 -192 -415 -237 -132 - 350 -161 -557 -210 
50) -13] -211 -1322 -355 -279 ~476 ~264 -332 =327 

0 30 -101 -642 -6 93 -26 52 -22 31 
25 -94 -170 <-1210 -226 -113 -255 -143 -692 -178 
50 -152 -211 -1477 -396 -257 -397 -277 -323 -340 
75 +-162 -234 -1643 -~522 -414 -582 -391 -412 -484 
100 -192 -262 -1839 -650 -598 -748 -518 -571 -638 
110 -197 -267 -1937 -697 -698 -827 -573 -~641 -713 
120 -208 -274 -2009 -721 -768 -880 -603 -692 -759 
130 -212 -277 -2061 -752 -839 -941 ~643 =751 -810 
140 -240 -281 -2173 -789 -974 -1005 -706 -848 -~878 
150 -251 -268 -2228 -785 -1052 -1026 -731 -912 -902 
160 -285 -256 -2338 -763 -1191 -1068 -769 -1026 -907 
170 -368 -191 -2508 -635 -1476 -1077 -791 -1280 -794 
175 -418 -112 -2589 -470 -1708 -1013 -774 -1506 -602 
180 -686 2218-276) 100 -2483 -671 -704 -2285 1062 
183 -2825 
2. Intercept at P=0 using values at 25 and 50 psi. 

At Be Cc D E F G H J 
INT -36 -129 -943 -56 31 -113 -9 -61 -16 
(-192@25) 

3. Zeroed Circumferential Strain Data (using intercept). 
Pres Ar Be Cc D E F G H J 

0 0 0 0 0 0 0 0 0 0 
Z5 -56 -63 528 -181 -163 -237 -152 -496 -194 
50 -95 -82 -379 -299 -310 - 363 -255 -271 -311 

0 0 0 0 0 0 0 0 0 0 
25 -58 -41 -267 -170 -144 -142 -134 -131 -162 
50 -116 -82 -534 -340 -288 -284 -268 -262 -324 
75 -126 -105 -700 - 466 -445 -469 - 382 -351 -468 
100 -156 -133 -896 -594 -629 -635 -509 -510 -622 
110 -161 -138 -994 -641 -729 -714 -564 -580 -697 
120 -172 -145 -1066 -665 -799 -767 -594 -631 -743 
130 -176 -148 -1118 -696 -870 -828 -634 -690 -794 
140 -204 -152 -1230 -733 -1005 . -892 -697 -787 -862 
150 -215 -139 .-1285 -729 -1083 -913 -722 -851 -886 
160 -249 -127 -1395 -707 -1222 -955 -760 -965 -891 
170 +332 -62 -1565 -579 -1507 -964 -782 -1219 -778 
175 -382 17 -1646 -414 -1739 -900 -765 -1445 -~586 
180 -650 356 -1818 156 -2514 -558 -695 -2224 1078 
183 -2856 


135 


K 
=42 
=2206 
- 367 
-13 
=20)3 
-413 
-448 
-543 
-600 
-629 
=0605 
-726 
=782 
-839 
-~976 
-1104 
-~ 1606 


-233 
-~374 


=210 
-420 
-455 
-550 
-607 
-~636 
-672 
=133 
-789 
-~846 
-983 
=1111 
=)o14 


—= 


i cE te oi ee 





STRAIN GAGE 


DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 


cylinder # 


18 


layup 
(902/(+_30)2]s 


1. Raw Circumferential Strain Gage Data. 


Pres 
0 
25 
50 
75 
100 
110 
20 
130 
140 
150 
160 
170 
175 


As 
-10 
a1e32) 
-82 
= 
= 1 
ae 7 
-203 
-24]1 
23 
- 366 
-485 
-787 


Be 
-6 
~42 
-65 
-83 
~94 
= 2l7) 
-~78 
-64 
-46 
-15 
71 
263 


2. Intercept at 


INT 


Asx 
-40 


Bsa 


-19 


c 
4 
-26 
-112 
=u Sh) 
-427 
-509 
-602 
~709 
-841 
= 99 
-1249 
-1815 


P=0 using values at 25 and 50 


c 
60 


D 
2 
-106 
-1690 
-186 
=-179 
-145 
-104 
-32 
69 
229 
545 
1314 


D 
-52 


E 

S| 
-285 
-580 
=914 
-1298 
-1485 
-1691 
-1929 
-2221 
-2576 
~3111 
~4174 
-5200 


5 
10 


F 
2 
79 
1260 
1100 
1000 
993 
975 
969 
967 
976 
1009 
1067 


F 
-1102 


3. Zeroed Circumferential Strain Data (using intercept). 


Pres 
0 
Z> 
50 
75 
100 
110 
120 
130 
140 
10) 
160 
170 
175 


Ar 
0 
-21 
~42 
-77 
-121 
-137 
-163 
-201 
-253 
- 326 
-445 
-747 


Be 
0 
-23 
-46 
~64 
-75 
-68 
-~59 
-45 
-27 


90 
282 


C 
0 
-86 
-172 
-319 
-487 
=209 
-662 
~769 
-901 
Os 
=1309 
-1875 


D 
0 
~54 
~108 
-~134 
-127 
=—33 
-52 
20 
121 
281 
597 
1366 


E 
0 
~295 
-590 
-924 
~ 1308 
-1495 
-1701 
-1939 
=2231 
- 2586 
~3121 
-4184 
~5210 


[2G 


F 

0 
1181 
2362 
2202 
2102 
2095 
2077 
2071 
2069 
2078 
2111 
2169 


test date 
12/11/91 

G H 

-3 -~8 
-153 -155 
-306 - 306 
-447 -448 
-588 -609 
-~634 -~678 
-~677 -762 
-711 -867 
-733 -1011 
-724 -1199 
-631 -154) 
-269 -2385 
psi. 
G H 

0 ~4 
G H 

0 0 
-153 -151 
~306 -302 
-447 -444 
-588 -605 
-634 -674 
~677 -758 
-711 -863 
-733 -1007 
~724 -1195 
-631 =+-1537 
-269 -2381 


J 


-143 
-259 
~348 
~425 
-433 
-445 
-439 
-416 
— Se) 
=191 
280 


= 7h 


-116 
~232 
~321 
-398 
-406 
-418 
-412 
=309 
- 324 
-164 
307 


K 


-164 
-296 
-414 
-536 
-566 
=612 
-661 
=721 
-784 
~869 
ees 


-32 


-132 
~264 
- 382 
-504 
-534 
-580 
~629 
-~689 
-752 
~ 837 
~944 





STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 


cylinder # layup test date 
23 (902/0)4t 11/22/91 


1. Raw Circumferential Strain Gage Data. 


Pres Asx Be c D E F G H J K 

0 Sep 331 10 -4 28 26 11 -9 5 -6 
25 -96 -414 -82 -92 -96 107 -99 -89 -134 -76 
50 -147 -471 -88 -218 -175 81 -163 -119 -241 -112 
75 -191 4-529 -68 - 321 -~254 5 -245 -144 -387 -142 
100 -199 -581 -104 -436 -303 -77 -312 -137 -601 -139 
110 -199 -601 -126 -480 -299 -97 -329 -106 -686 -99 
120 -211 4-635 -171 -565 -273 -143 -346 -55 -805 - 33 
130 -188 -661 -203 -645 -209 -233 -363 25 -1037 76 
140 -131 -729 -332 -857 50 -488 -350 263 -1524 430 
150 59 -926 -834 -1306 1028 -1305 -107 849 -2559 1362 
2 | -4400 


2. Intercept at P=0 using values at 25 and 50 psi. 


Ax Be Cc D E F G H J K 
INT -45 -357 -76 34 -17 133 -35 -59 -27 -~40 


3. Zeroed Circumferential Strain Data (using intercept). 


Pres Ar Be c D E F G H J K 

0 0 0 0 0 0 0 0 0 0 0 
25 -51 -57 -6 -126 -79 -26 -64 -30 -107 -36 
50 -102 -114 -12 -252 -158 -52 -128 -60 -214 -72 
75 -146 -172 8 - 355 -237 -128 -210 -85 - 360 -102 
100 -154 -224 -28 -470 -286 -210 -277 -78 -574 -99 
110 -154 -244 -50 -514 -282 -230 -294 -47 -659 -59 
120 -166 -278 -95 -599 -256 -276 -31] 4 -778 i 
130 -143 -304 -127 -679 - 192 -~ 366 -328 84 -1010 116 
140 -86 -372 -256 -891 67 -621 -315 322 -1497 470 
150 104 4-569 -758 -1340 1045 -1438 -72 908 -2532 1402 
152 - 4373 


i Ry LE ot oe 


APPENDIX E 


ABD STIFFNESS MATRICES 


138 





of t. 


((90)4/(0)2)6) r = 2.755", 1 = 16" 


Buckling Pressure for a cylinder with dinensions of: 


r= 2.75, 


Per = 
ABAQUS Pcr = 


ABD Matrix: 
§.3129E+05 


3.0428E+04 
0.0000 


0.0000 
0.0000 


1 = 18.00, t = 0.0720 => Zs 
and a load ratio of; Ny = 2.00 Nx 
218.2 NX = 300.6 
242.5 
3.0428E+04 0.0000 0.0000 
9.5390E+05 0.0000 0.0000 
0.0000 7.2000E+04 0.0000 
a eee | co 
0.0000 0.0000 67.23 
0.0000 0.0000 13.15 
0.0000 0.0000 0.0000 








((90)3/(0)3)s; r = 2.755", 1 = 186" 


1633. 
= 1.00 


M = 1 N 


0.0000 
0.0000 
0.0000 


0.0000 


Buckling Pressure for a cylinder with dimensions of: 


r= 2.75, 


Pcr = 
ABAQUS Pcr = 


ABD Matrix: 
7.4260E+05 


3.0428Et04 
0.0000 


0.0000 
0.0000 





1 = 18.00, t = 0.0720 => Zs 
and a load ratio of: Ny = 2.00 Nx 
236.6 Nx = 325.9 
230.5 
3.0428E+04 0.0000 0.0000 
7.4260E+05 0.0000 0.0000 
0.0000 7.2000EB+04 0.0000 
"7 ses 2 2e« se» es& 2 @ -= ¢ 2"#®= = = @ 
0.0000 0.0000 115.4 
0.0000 0.0000 13.15 
0.0000 0.0000 0.0000 


0.0000 





((90)2/(£30)2)s; r = 2.755", 1 = 18" 


1633. 
= 1.00 


M = 1 N 


0.0000 
0.0000 
0.0000 
13.15 
526.2 
0.0000 


Buckling Pressure for a cylinder with dimensions of: 


r= 2.75, 1 18.00, ¢ = 0.0720 => ts 1633. 
and a load ratio of: Ny = 2.00 Nx = 1.00 
Per = 206.6 Nx = 284.6 M = 1 N 
ABAQUS Pcr = 205.4 
ABD Matrix: 
6.0055E+05 1.7247E+05 0.0000 0.0000 0.0000 
1.7247E+05 6.0055E+05 0.0000 0.9000 0.0000 
0.0000 0.0000 2.1404E+05 0.0000 0.0000 
- = 2 se ww 2@ @ w» @ @ =e =& 2 @B & & & @& +$ eeewrwreerr wf FrFeeer-e s-=f#ee 
0.0000 0.0000 0.0000 141.4 40.42 
0.0000 0.0000 0.0000 40.42 445.7 
0.0000 0.0000 0.0000 19.08 7.271 








, rf/t = 38. 


ad 2 


0.0000 


0.0000 
0.0000 


, r/t = 38. 


- 2 


0.0000 
0.0000 
0.0000 


r/ft = 38. 


= 3 


0.0000 


0.0000 
0.0000 


- 
Se 


ee 


a 





Buckling Pressure for a cylinder with dimensions of: 
rs ao, 1 «= 18.00, t = 0.0720 => Z = 1633, r) r/t = 38. 
and a load ratio of: Ny = 2.00 Nx = 1.00 


Pcr = 184.1 Nx = 253.7 M = 1 N= 3 
ABAQUS Pcr = 178.8 


ABD Matrix: 





5.3129E+05 3.0428E+04 0.0000 2536. 0.0000 0.0000 
3.0428E+04 9.5390E+05 0.0000 0.0000 -2536. 0.0000 
0.0000 0.0000 7.2000E+04 0.0000 0.0000 0.0000 
~ eee eee eee eee eee ~-¢2e 22226222 22 © e ee ew ew ew 
2536. 0.0000 0.0000 234.6 13.15 0.0000 
0.0000 -2536. 0.0000 13.15 407.0 0.0000 
0.0000 0.0000 0.0000 0.0000 0.0000 31.10 





Jones (90,90,90,60,30,0,0,15,30,60,75,90)? r = 2.755", 1 = 18" 
Buckling Pressure for a cylinder with dimensions of: 
r= 2.75, 118.00, t= 0.0720 => Z= 1633. , v/t = 38, 
and a load ratio of: Ny = 2.00 Nx = 1.00 


Pcr = 235.9 Nx = 324.9 M = 1 N= 3 
ABAQUS Pcr = 213.4 


ABD Matrix: 





5.5408E2+05 1.1329Z+05 1.1791E+05 1286. 426.1 565.4 
1.1329E+05 7.6539E+05 1.1791E+05 426.1 -2138. 934.5 
1.1791E+05 1.1791E+05 1.5486E+05 565.4 934.5 426.1 
ate a a = as = = eee ee $e ee we ew we we we we = ow - «= « 
1286. 426.1 565.4 83.12 35.45 22.95 
426.1 -2138. 934.5 35.45 513.9 43.62 
565.4 934.5 426.1 22.95 43.62 53.40 





nth 


‘ 


ee Se wee ee eee 
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